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FOREWORD 


In  June  1957,  Armour  Research  Foundation  completed  a  feasi¬ 
bility  study  of  a  lightweight  field  artillery  weapon  intended  for  .use  by  the 
Marine  Corps.  The  study  specifically  explored  the  possibility  of  develop¬ 
ment  of  a  helicopter-transportable  field  weapon  with  a  high- rate -of-fire 
burst  capability.  The  results  of  this  study  indicated  that  a  lightweight  field 
artillery  weapon  could  be  developed  to  use  boosted  rocket  ammunition,  be 
capable  of  helicopter  transport,  and  fire  six  rounds  in  approximately  two 
seconds.  A  project  was  then  initiated  at  ARF  to  design  and  develop  such  a 
launcher. 

The  first  phase  of  the  program  consisted  of  the  design  and 
fabrication  of  the  XM70  Launcher,  Prototype  No.  1.  This  prototype  dis¬ 
played  an  exceptionally  optimistic  performance  such  that  a  second  phase  of 
the  program  was  initiated.  The  second  phase  consisted  of  the  incorporation 
design  improvements  and  construction  of  two  more  prototypes,  XM70E1, 
Prototypes  2  and  3.  Subsequently,  Prototypes  4  and  5,  XM70E2,  were 
designed  and  built. 

Although  the  early  performance  of  the  XM70  Launcher  was  satis¬ 
factory,  it  was  realized  that  improvements  in  performance  were  possible 
and  necessary.  In  order  to  determine  the  best  methods  of  improving  perform 
ance,  specifically  the  launcher  dynamic  response  and  its  effect  upon  accuracy 
a  long-range  study  program  >vas  initiated.  This  program  consisted  of  a 
concentrated  effort  of  both  experimental  and  analytical  investigations  of  the 
launcher  dynamics. 
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This  long-range  study  program  was  begun  in  October  of  1959  and 
concluded  in  March,  1961,  Progress  reports  concerning  this  program  were 
included  in  the  bi-monthly  reports  of  the  design  and  development  project. 
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ABSTRACT 

Prototype  No.  3  Launcher,  XM70E1,  115  mm  was  instrumented  wi 
strain,  pressure  and  displacement  gages;  these  furnished  the  actual  loading 
and  motion  of  the  launcher  structure.  In  addition  to  certain  simple  dynamic  • 
analyses,  a  3-dogree-of-freedom,  nonlinear  mathematical  model  of  the 
launcher  dynamics  was  derived  and  programmed  for  solution  on  Armour 
Research  Foundation's  Univac  1105.  The  output  of  the  computer  program 
was  correlated  with  experiment  and  used  to  study  the  effect  of  physical 
parameter  variations. 

Regions  of  instability  of  the  launcher  motions  were  shown  to 
exist  for  burst  firings;  relationships  between  component  stiffness  and 
damping  were  found  which  optimized  the  launcher  response  to  firing  loads, 
based  upon  a  simple  accuracy  criterion.  Certain  design  suggestions  were 
evaluated  and  shown  to  benefit  the  accuracy  of  the  launcher. 
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LONG-RANGE  STUDY  PROGRAM 
LIGHTWEIGHT  ARTILLERY  WEAPON 


I.  INTRODUCTION 

A  unique  automatic  rocket  launcher  in  the  field  artillery  weapon 
class  has  been  designed  and  developed  by  Armour  Research  Foundation 
under  the  technical  supervision  of  Rock  Island  Arsenal.  This  weapon  is 
a  115mm  launcher  capable  of  firing  single  shots  or  bursts  of  up  to  six  rounds. 
Boost  charge  variations  combined  with  rocket  thrust  increments  and  variable 
elevation  provide  an  exceptionally  wide  range  capability.  The  launcher  can 
be  traversed  +,  20°  without  disturbing  its  original  emplacement.  Based  upon 
early  feasibility  studies  and  also  upon  the  performance  of  the  first  proto¬ 
type,  this  lightweight,  automatic  launcher  concept  showed  itself  to  be 
unusually  capable  as  an  effective  weapon.  It  was  realized,  however,  that 
production  schedule  requirements  for  succeeding  prototypes  might  hinder 
the  development  of  the  concept,  because  of  the  tendency  to  freeze  important 
design  features.  As  a  result,  this  long-range  study  program  was  initiated 
to  investigate  the  dynamic  behavior  of  the  launcher  during  burst  firing, 
relate  this  behavior  to  accuracy,  and  determine  the  important  physical 
parameter  values  which  will  insure  good  accuracy. 

The  program  was  comprised  of  theoretical  and  experimental 
studies  of  the  dynamics  of  the  launcher.  These  consisted  of  finding  the 
effects  of  certain  parameter  variations  upon  the  launcher  motion;  accuracy 
was  related  to  the  motion  by  simple  criteria.  Some  of  the  important  launcher 
parameters  which  were  considered  are:  structural  rigidity,  damping, 
recoil  system  characteristics,  launcher  geometry,  and  ground  restraints. 

The  majority  of  the  work  accomplished,  both  theoretical  and 
experimental,  is  presented  in  the  appendices.  A  description  of  the  work  and 
results  as  well  as  the  important  conclusions  are  presented  in  the  following 
sections.  Section  II  is  an  introductory  portion  of  the  report  containing  a 
description  of  the  launcher  structure  and  its  relationship  to  its  dynamic 
response.  Some  of  the  operational  characteristics  of  the  launcher  system, 
as  they  affect  the  response,  are  also  discussed. 
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Section  III  presents  the  results  of  the  experimental  work,  and  is 
broken  into  three  parts.  These  respectively  describe  stiffness  measurements 
of  the  structural  components  and  two  separate  firing  programs.  Stiffness 
measurements  were  made,  both  statically,  of  all  the  major  structural  com¬ 
ponents,  and  dynamically  of  the  entire  supporting  structure  assembly  (carriages 
and  trails).  The  second  experimental  effort  was  a  firing  program  wherein 
approximately  75  proof  slugs  were  fired.  During  these  firings  the  launcher 
was  instrumented  with  pressure  gages,  displacement  transducers,  and  strain 
gages  to  provide  information  on  the  force  input,  the  response  of  the  launcher, 
and  accuracy  in  the  vertical  plane.  Finally,  horizontal  dispersion  was 
significantly  reduced,  as  shown  by  targets,  with  the  use  of  two  fixtures  -- 
one  restraining  the  trails  from  horizontal  motion  while  the  second  stiffened 
the  recoiling  assembly  framework. 

In  Section  IV,  the  theoretical  analyses  are  described;  they  consist 
of  the  following  mathematical  models: 

1.  A  planar,  nonlinear  3-degree-of-freedom  model  (called 
Mathematical  Model  V  in  progress  reports), 

2.  A  third-order,  linear  model , 

3.  A  linear,  2-degree-of-freedom  model , 

4.  A  three-dimensional .  nonlinear,  n-degree-of-freedom  system, 

5.  A  second-order,  linear  system. 

The  first  of  these  was  implemented  with  a  computer  program  and  was  widely 
applied  to  study  launcher  motion.  It  contained  representations  of  the  hydro¬ 
pneumatic  recoil  system,  the  indexing  cams,  and  flexibilities  and  damping 
of  both  the  supporting  structure  and  elevating  system.  Information  about 
the  response  of  linear,  third-order  systems  led  to  a  structural  optimization 
of  the  elevating  system.  The  2-degree-of-freedom  system  was  an  approxi¬ 
mation  of  the  launcher  used  to  determine  the  best  location  for  auxiliary 
dampers.  The  fourth  model  is  a  relatively  complex  representation  of  the 
launcher;  its  significant  contribution  was  the  inclusion  of  structural 
flexibilities  of  both  the  recoiling  assembly  framework  and  firing  tube  and  also 
consideration  of  motion  out  of  the  vertical  plane.  Unfortunately,  this  model 
could  only  be  developed  to  the  point  where  it  was  programmed  for  the  Univac 
but  not  run.  The  final  system  was  a  single-degree-of-freedom  system  used 

unoun  :  S  E  S  E  A  B  C  W  P  O  U  N  0  *  T  I  O  N  OP  ILLINOIS  INSTITUTE  OP  TECHNOLOGY 


to  investigate  the  important  parameter  values  which  can  cause  the  response 
to  decay  by  a  desired  amount  between  rounds  of  a  burst. 

The  3-degree -of-freedom  model  was  used  to  correlate  theoretical 
response  curves  with  the  experimental  response  obtained  with  various  structural 
changes.  The  correlation  was  very  good  for  certain  parameter  combinations, 
but  the  theoretical  response  curves  lacked  certain  characteristics  as  some 
parameters  were  varied.  The  reason  for  this  is  postulated;  allowing 
conclusions  to  be  drawn  which  indicate  the  desirability  of  certain  combinations 
of  structural  parameters  which  correspond  to  improved  accuracy. 

A  theoretical  investigation  was  made  of  a  specific  design  recommenda¬ 
tion  --  delayed  recoil  rods.  Because  the  3 -degree -of-freedom  model  had 
an  explicit  representation  of  the  recoil  Bystem,  the  behavior  resulting  from 
increasing  the  orifice  area  during  the  time  the  rocket  is  in  the  launcher 
could  be  studied.  This  cutback  on  the  recoil  rods  resulted  in  near  free  recoil 
for  a  short  time. 

The  final  item  studied  was  done  as  a  matter  of  interest;  the  terms 
representing  the  Coriolis  acceleration  of  the  recoiling  assembly  were 
neglected  during  one  computer  run.  Comparison  was  then  made  with  a 
normal  response  to  determine  the  effect. 

Although  the  launcher  has  a  +  20°  traverse  capability,  all  of  the 
studies  considered  only  zero  traverse  accuracy  because  of  its  order  of  impor¬ 
tance.  Ground  restraints  and  soil  considerations  were  included  in  the  n-degree- 
of-freedom  model  and  could  not  be  studied  because  of  the  status  of  the  model. 
One  method  of  auxiliary  ground  restraint,  horizontal  trail  ties,  was  tried  and 
is  discussed  with  horizontal  dispersion  reduction. 
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II.  DESCRIPTION  OF  THE  LAUNCHER  STRUCTURE 

The  complete  launcher  is  shown  in  two  firing  positions,  in 
I  and  2.  It  consists  of  two  main  assemblies:  the  tipping  assembly  and 
the  supporting  structure.  The  total  weight  of  the  launcher,  3150  lb,  is  almost 
equally  divided  between  these  two  assemblies.  The  tipping  assembly  rotates 
relative  to  the  supporting  structure  at  the  trunnion;  the  elevating  system 
separates  them  at  any  fixed  position-from  -6°  to  65°. 

A.  Tipping  Assembly 

The  tipping  assembly  is  composed  of  a  cradle  and  recoiling 
assembly;  the  former  is  simply  the  containment  for  the  latter.  The  firing 
tube  is  fixed  to  the  front  of  the  structural  framework  of  the  recoiling  assembly 
(see  Fig.  3);  it  also  supports  two,  contrarotating  breech  clusters  of  three 
rounds  each  and  holds  the  recuperator  and  recoil  system.  The  framework 
consists  essentially  of  the  front  plate,  rear  plate,  two  side  tie  bars,  the 
recuperator  cylinder  and  a  lower  tie  bar;  the  recoil  cylinders  are  not  load 
bearing  members  of  this  framework.  Two  indexing  cams  extend  into  the 
hubs  of  the  breech  clusters  and  cause  them  to  alternately  rotate  into  firing 
position  during  counterrecoil.  These  cams  are  fixed  to  the  front  of  the 
cradle. 

B.  Supporting  Structure 

The  supporting  structure  is  the  name  given  to  the  portion  of  the 
launcher  below  the  trunnions;  this  structure  is  made  up  of  what  are  commonly 
called  the  upper  carriage,  lower  carriage  and  trails.  An  exploded  view  of 
the  supporting  structure  in  Fig.  4  shows  its  essential  components.  This 
entire  structure  is  composed  of  hollow,  sheet  steel  construction  with  various 
necessary  non-structural  components  accounting  for  20%  of  its  total  weight. 
These  components  are  the  traverse  and  eievating  mechanisms,  wheels,  firing 
base,  lunnette  and  a  jack  for  raising  and  lowering  the  firing  base.  Ir.  firing 
position,  the  trails  are  separated  by  an  included  angle  of  40°  which  permits 
+  20°  traverse  firing;  the  upper  box,  trunnion  sides  and  tipping  assembly 
rotate  relative  to  the  remaining  structure  to  attain  the  desired  traverse 
positions. 
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C.  Elevating  System 

The  elevating  system  plays  an  important  role  in  the  dynamic 
response  of  the  launcher  and  deserves  some  discussion.  This  system  is 
composed  of  hollow  threaded  rods  passing  through  ball  bearing  nut  assemblies 
held  in  brackets  on  each  side  of  the  cradle.  The  lower  ends  of  these  rods  are 
inserted  into  the  trunnion  sides  where  they  are  supported  by  a  thrust  bearing 
in  a  gear  box.  The  rods  themselves  are  relatively  stiff  in  tension  and 
compression;  certain  other  components  of  this  system,  however,  comprise 
a  significant  flexibility.  The  brackets  at  the  cradle  are  cantilevered  and 
are  inherently  flexible;  torsion  of  the  cradle  sides  also  contributes  to 
effective  relative  motion  of  the  tipping  assembly  and  supporting  structure. 

The  ball  bearings  in  the  screw-nut  mechanism,  due  to  unavoidable  clear¬ 
ances  also  add  to  the  system  flexibility  in  the  form  of  nonlinearities.  Like¬ 
wise,  the  lower  connections,  through  local  deflections  of  the  sheet  metal  box 
and  trunnion  sides,  are  relatively  flexible.  Thus,  significant  rotation  of  the 
tipping  assembly  mass  center  can  occur  about  the  trunnions  (in  the  vertical 
plane)  during  firing, due  to  the  flexibility  of  the  elevating  system. 

D.  Damping 

As  designed,  damping  of  the  launcher  motions  can  occur  from 
four  sources: 

1.  Structural  material  damping, 

2.  Interaction  with  the  soil, 

3.  The  hydraulic  recoil  system, 

4.  The  support  buffers. 

The  first  damping  source  is  negligible  in  the  presence  of  the  others.  The 
second  source  varies  with  the  type  of  soil  and  can  be  either  negligible  or 
significant.  The  third  source  of  damping  is  important  because  relative  motion 
between  the  cradle  and  recoiling  assembly  can  cause  relatively  large  amounts 
of  energy  to  be  damped  out  by  forcing  oil  through  the  orifice  in  each  recoil 
cylinder.  The  support  buffers  consist  of  two  cylinder  and  piston  assemblies 
attached  to  the  front  of  the  lower  carriage;  the  cylinders  are  fixed  to  the 
carriage  and  the  piston  rods  rest  against  the  firing  base.  When  the  supporting 
structure  rotates  forward,  about  the  ball  joint  connection  to  the  base,  the 
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pistons  force  oil  through  orifices  and  damp  out  this  forward  motion.  Small 
springs  in  parallel  with  the  pistons  restore  the  position  of  the  supporting 
structure.  These  buffers  are  mainly  to  prevent  the  launcher  from  tipping 
over  following  the  final  shot  of  a  burst  when  the  recoiling  assembly  returns 
to  battery.  These  buffers  operate  during  bursts,  however,  when  any  forward 
rotation  of  the  supporting  structure  occurs. 

E.  Deflection 

Before  discussing  the  dynamic  response  of  the  launcher  structure 
to  firing  loads,  it  is  necessary  to  first  examine  the  supporting  structure  to 

see  how  it  deflects  to  a  load  at  the  trunnions.  We  will  consider  zero  traverse, 

•  * 

a  completely  symmetric  structure  and  loads  in  the  vertical  plane  only.  Let 
us  assume  that  the  tipping  assembly  and  elevating  system  are  removed;  a 
force  at  the  trunnions  (for  example,  at  45°  with  respect  to  the  ground)  causes 
the  trunnion  sides  to  deflect  relative  to  the  upper  box.  Since  these  members 
are  essentially  short,  stubby  beams,  both  bending  and  shear  are  of  equal 
significance.  We  now  have  a  shear  force  and  moment  at  the  sides  of  the 
upper  box;  this  causes  bending  and  shear  of  both  sides  of  this  member  rela¬ 
tive  to  the  bearing.  However,  th&  loading  also  causes  torsion  of  the  upper 
box  section,  and  is  the  most  significant  deflection  of  this  member  compared 
to  shear  and  bending  deflections.  The  shear  load  is  further  transferred 
through  the  bearing  down  to  the  ground  at  the  firing  base.  The  moment  reaction, 
however,  is  transferred  through  the  bearing  and  causes  torsion  of  the  lower 
box,  which  eventually  places  a  bending  moment  on  the  front  of  each  trail.  With 
respect  to  trunnion  deflection  in  the  vertical  plane  we  then  see  that  the 
following  values  of  stiffness  of  the  structure  components  are  important: 

1.  Bending  and  shear  of  the  trunnion  sides, 

2.  Torsion  of  the  upper  and  lower  box, 

3.  Separation  of  the  traverse  bearing, 

4.  Bending  of  the  trails. 

Significant  deflections  occur  both  from  the  supporting  structure  and  from  the 
elevating  system  and  cause  a  nonlinear  relationship  between  the  various 
deflections, due  to  the  geometry. 
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F •  Experimental  Structure  Modifications 

During  the  experimental  firings  conducted  on  this  program, 
various  modifications  to  the  structure  were  made.  These  took  essentially 
two  forms,  changes  in  stiffness  and  changes  in  damping.  Modifications  of 
the  structure  which  essentially  affect  only  the  response  in  the  vertical  plane 
are  discussed  first.  The  stiffening  fixtures  are: 

1.  Variable  stiffness  elevating  rods, 

2.  Upper  and  lower  box  end-clamps, 

3.  Trunnion-to -trail  struts. 

The  utility  of  the  first  fixture  is  self  explanatory.  The  second  and  third 
sets  of  fixtures  were  designed  to  permit  a  temporary  increase  in  the  stiff¬ 
ness  of  the  supporting  structure  by  decreasing  the  torsion  of  the  upper  and 
lower  boxes.  The  clamps  were  1/2  in.  thick  steel  brackets  which  joined 
each  of  the  four  corners  of  the  boxes  and  transferred  the  moment  reaction 
from  the  base  of  the  trunnion  sides  directly  to  the  trails.  The  trunnion-to- 
trail  struts  caused  most  of  the  trunnion  reaction  to  be  directly  transferred 
to  the  trails.  For  these  struts,  the  trails  had  to  be  rotated  inward  until  they 
were  set  directly  beneath  the  trunnions.  Having  the  trails  straight  back 
caused  an  additional  stiffness  increase  in  this  configuration. 

Damping  was  added  to  the  elevating  system  in  either  of  two  ways. 
Auxiliary,  hydraulic  dampers  were  placed  between  the  cradle  front  and 
trunnion  sides  in  parallel  to  the  elevating  screws.  Figure  5  shows  these 
dampers  on  the  launcher.  Provisions  were  made  to  change  the  damping  source 
by  the  iMe  of  frictional  ring  springs  in  the  variable  flexibility  elevating  rods. 

A  more  complete  description  of  these  fixtures  is  presented  in  Appendix  G. 

The  above  fixtures  were  intended  to  cause  a  direct  change  of  the 
launcher  response  in  the  vertical  plane  (neglecting  coupling  of  vertical  and 
horizontal  motions).  Another  temporary  change  to  the  launcher  structure 
was  made,  this  time  for  the  purpose  of  affecting  the  horizontal  motion  of  the 
launcher.  This  change  was  the  addition  of  a  thin,  sheet  metal  shell  which 
surrounded  the  recoiling  assembly  and  was  attached  to  the  side  tie  bars. 

This  shell  stiffened  the  framework  in  a  manner  such  that  the  front  plate  and 
rear  plate  would  be  more  securely  restrained  from  rotating  with  respect  to 
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each  other  about  a  vertical  axis.  This  type  of  rotation  was  experienced  in 
the  launcher  and  caused  a  significant  horizontal  dispersion.  For  the  same 
purpose,  reduction  of  horizontal  dispersion,  the  trail  ends  were  fastened 
by  cables  to  small  stakes  at  the  side  of  each  trail.  This  prevented  the  trails 
from  rotating  (in  the  plane  of  the  ground)  and  caused  a  decrease  in  horizontal 
dispersion. 
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III.  EXPERIMENTAL  WORK 

The  experimental  effort  was  composed  of  three  separate  ventures: 

1.  The  determination  of  the  structural  component  flexibilities, 

2.  Determination  of  the  dynamic  response  of  the  launcher  in  the 
vertical  plane, 

3.  Testing  of  fixtures  whicli  were:to  decrease  horizontal  dispersion. 
The  values  of  stiffness  of  the  structural  components  were  needed  for  two 
reasons:  first,  as  input  to  the  theoretical  studies  and  second,  to  determine 
if  the  particular  structural  designs  yielded  the  desired  stiffness,  as  required 
by  earlier  analyses.  The  stiffness  was  found  experimentally  because  of  the 
complexity  of  the  structure;  it  was  more  ■effieient  than>corresporidix\g  \ 
calculations.  In  das'4st whirl  calculations  were  made,  the  measurements 
afforded  a  method  of  checking. 

The  dynamic  response  in  the  vertical  plane  was  determined  for 
correlation  with  theory  and  thereby  used  to  validate  the  corresponding 
mathematical  model,  and/or  point  out  its  deficiencies. 

A.  Physical  Parameter  Measurements 

Static  stiffness  measurements  of  the  following  components  or 
assemblies  were  made: 

1.  The  entire  supporting  structure, 

2.  The  entire  supporting  structure  with  the  upper  and  lower  boxes 
clamped, 

3.  The  entire  supporting  structure  with  the  box  clamps  and  with  the 
trunnion-to-trail  struts. 

4.  The  entire  supporting  structure  in  traverse  position, 

a.  20°  right 

b.  20°  left 

5.  The  lower  elevating  screw  support, 

6.  Trails, 

7.  Elevating  screw  -  upper  bracket  combination. 

Since  rotational  motion  is  the  sole  significant  contributor  to 
dispersion,  the  moment  of  the  recoil  force  per  unit  rotation  of  a  line  through 
the  ball  joint  and  the  trunnions  (in  the  vertical  plane),  is  used  as  an 
indication  of  supporting  structure  stiffness.  The  static  stiffness  of  the 
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supporting  structure  is  therefore  given  in  units  of  in.  -lb/ radian.  For  the 
entire  supporting  structure  at  zero  traverse,  this  was  measured  to  be 
32  x  10^  in.  -lb/ radian.  When  the  box  clamps  were  added,  this  stiffness 
became  38.  5  x  10^,  approximately  20%  greater.  A  much  larger  increase 
occurred  when  the  trunnion-to -trail  struts  were  added;  the  stiffness  rose  to 
55,  5  x  10^,  a  73%  increase.  This  was  due  not  only  to  the  direct  load  transfer 
from  the  trunnion  to  the  trails  but  is  also  because  the  trails  were  straight 
back;  in  this  position  they  contributed  more  rigidity  to  moments  in  the  vertical 
plane.  During  traverse  (20°  left  or  right)  tlfe  stiffness  of  the  supporting 
structure  dropped  to  25  x  10^  in.  -lb/ radian  for  the  "as-designed"  structure. 
This  decrease  is  expected  since  the  stiffness  of  one  trail  is  almost  completely 
lost  in  this  configuration. 

The  trail  stiffness  (total,  for  both  trails  straight  back)  for  a 
bending  moment  applied  at  the  forward  end  resulted  in  145  x  10^  in.  -lb/ radian. 
In  this  case,  the  stiffness  was  lower  than  predicted  by  calculation,  and  is 
believed  due  to  local  deflections  in  the  area  of  the  trail  hinge  pins. 

Measurement  of  the  stiffness  of  the  lower  elevating  screw  supports 

illustrated  an  unanticipated  result.  The  left  support  was  nonlinear  up  to  a 

load  of  800  lb;  at  this  load  it  assumed  a  relatively  linear,  load-deflection 

5 

curve  shape,  with  a  slope  of  3.  5  x  10  lb/in.  The  right  support  had  a  linear 

5  * 

stiffness  of  3.  5  x  10  lb/in.  from  zero  load.  This  would  indicate  that 
unsymmetrical  elevating  system  deflections  would  occur.  This  actually 
happened  and  can  be  seen  in  the  experimental  records  in  Appendix  J. 

The  elevating  rods  with  their  upper  brackets  attached  exhibited 
nonlinear  load -deflection  curves  as  assemblies.  In  this  situation,  the  non- 
linearities  are  due  to  either  or  both  of  the  following  reasons: 

1.  Clearances  between  the  ball-nut  and  the  ball  bearings. 

2.  Nonlinear  local  deflections  of  the  balls  and  ball-nut. 

The  clearances  can  cause  nonlinearities  because  as  more  balls  come  into  con¬ 
tact  at  higher  loads,  the  stiffness  increases  until  all  the  bails  are  in  contact. 
The  resulting  stiffness  curves  are  shown  in  Fig.  H- 12. 

In  addition  to  the  static  stiffness  measurements,  dynamic 
measurements  were  made  of  the  supporting  structure.  As  opposed  to  the 
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static  measurements  where  the  deflection  is  measured  under  a  very  slowly 
varying  load,  the  dynamic  stiffness  is  a  measurement  of  the  rms  value  of  a 
sinusoidally  varying  force  and  the  rms  value  of  the  displacement.  (Maximum 
values  also  could  have  been  used.)  This,  of  course,  varies  with  frequency. 

A  standard  600  lb  electrodynamic  vibration  shaker  was  used  to 
excite  the  launcher  structure  for  the  dynamic  stiffness  measurements. 

Figure  H~2c  shows  the  physical  arrangement.  The  tipping  assembly  and 
elevating  rods  were  removed  and  a  relatively  stiff  bar  was  placed  across  the 
trunnions.  The  shaker  was  coupled  to  the  bar  with  a  dynamic  force  trans¬ 
ducer  which  yielded  the  input  force.  Accelerometers  were  placed  on  the 
trunnions  to  yield  the  motion  ^output;  the  dymamic  stiffness,  k^,  is  simply, 


where 

F  =  Measured,  input  force,  rms  (or  peak), 

C7  =  Input  force  frequency, 

a  =  Measured,  linear  trunnion  acceleration,  rms  (or  peak). 

Although  the  repeatability  of  the  magnitude  of  the  measurements  was  not 
good,  the  shape  of  the  stiffness -frequency  curve  was  the  same  for  each  of  the 
two  runs.  The  relative  minimum  values  of  this  curve  (Fig.  H-ll)  correspond 
to  resonance  frequencies  of  the  structure.  The  first  mode  frequency  (the  fre¬ 
quency  at  the  first  resonant  point)  of  the  entire  supporting  structure  is  near 
70  cps.  Since  the  firing  frequency  of  a  burst  is  slightly  less  than  3  rounds  per 
second,  we  can  be  confident  that  an  undesirable  structural  resonance  will  not 
occur.  The  recoil  force  itself  could  possibly  contain  variations  at  a  frequency 
near  70  cps,  however,  previous  experimental  and  calculated  curves  show  that 
the  magnitude  of  the  recoil  force  does  not  have  these  high  frequency  variations. 
The  result  is  that  for  this  launcher,  calculation  of  the  first  mode  response 
should  be  sufficient  to  predict  the  response  accurately. 

A  rough  check  of  the  magnitude  of  this  frequency  can  be  made, 
combining  the  measured,  static  stiffness  with  the  result  of  a  calculation  of  an 
equivalent  structure  inertia.  The  equivalent  structure  inertia  was  found  by 
assuming  that  the  structure  and  its  equivalent  both  respond  harmonically 
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In  the  static  deflection  shape  of  the  actual  structure.  Imposing  the  conditions 
that  the  maximum  deflection  and  kinetic  energy  of  both  systems  are  identical 
yields  an  equivalent  inertia;  for  this  structure  it  is  115  in.  -lb-sec  ,  The 
first  mode  frequency  is  found  by  treating  the  equivalent  system  as  a  1 -degree - 
of-freedom  system.  This  frequency  is: 


f  =  84 

where 

k  =  measured  static  stiffness, 

8 

I  =  calculated  equivalent  inertia  about  the  ball  joint. 

These  two  frequencies  correspond  sufficiently  close  when  experimental 
errors  and  simplified  calculations  are  considered.  (The  methods  of  deter¬ 
mination  of  the  values  of  stiffness  of  the  launcher  structural  components  are 
described  more  fully  in  Appendix  H. ) 

B.  Launcher  Response  Determination 

The  dynamic  response  (both  displacement  and  velocity)  of  the 
launcher  in  the  vertical  plane  to  firing  loads  was  obtained  directly  by 
instrumenting  the  launcher  during  a  controlled  firing  program.  The  general 
procedure  was  to  fire  a  two  or  three  round  burst  for  each  structural  condition, 
followed  by  two  single  shots.  The  reason  six-round  bursts  were  not  necessary 
for  each  condition  is  that  the  response  to  rounds  3,  4.  5,  and  6  do  not  differ 
appreciably  from  the  second.  The  short  test  bursts,  therefore,  illustrated 
full  burst  conditions  and  the  single  shots  provided  a  reasonable  check  on 
repeatability  of  the  type  of  response.  The  firing  program  with  the  instrumen¬ 
tation  and  structural  configurations  listed  for  each  test  is  shown  in  Table  1-1. 
A  list  of  the  type  and  use  of  recorded  data  is  as  follows: 

1.  Strain  adjacent  to  the  firing  pin  (on  the  rear  plate  of  the  recoiling 
assembly) . for  a  firing  signal  proportional  to  the  boost  charge, 

2.  Strain  of  the  upper  trail  flanges  near  the  forward  pin  connections. to 
indicate  the  supporting  structure  rotation. 
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3.  Cradle  displacement,  in  the  vertical  plane,  at  the  tip  of  the  cradle. 
(Combined  with  similar  trunnion  motion,  this  furnishes  cradle 
rotation). 

4.  Elevating  rod  strain  to  furnish  the  force  in  the  rods, 

5.  Displacement  of  the  trunnion,  parallel  and  perpendicular  to  the 
cradle  (some  firings  only), 

6.  Damper-support  strain  to  furnish  the  damper  force, 

7.  Recoil  system  pressure  to  indicate  the  recoil  force. 

The  rotational  response  of  the  cradle  and  the  trail  strain  (structure 
rotation)  comprised  the  response  data  of  primary  significance.  The  cradle 
rotation  was  used  as  the  primary  indication  of  the  effectiveness  of  a  particu¬ 
lar  structural  change  in  reducing  dispersion.  The  trail  strain  formed  an 
important  quantity  in  correlating  theory  and  experiment.  The  forcing  functions 
in  the  form  of  rear  plate  strain  and  recoil  system  pressure  were  also  necessary 
quantities.  (It  should  be  mentioned  that  the  trail  strain  measured  is  linearly 
proportional  to  the  trail  rotation  only  during  positive  rotation,  the  initial  direc¬ 
tion  of  deflection.  This  is  because  during  negative  rotation,  the  supporting 
structure  rotates  forward,  the  trail  ends  leave  the  ground  and  the  strain  is 
then  a  measurement  of  the  trail  inertial  loads. ) 

Table  1  summarizes  the  launcher  conditions  during  the  various 
bursts  which  are  discussed  here.  Except  for  the  modifications  listed,  the 
launcher  was  as  designed. 

Table  1 


Burst 

No. 

Modifications  to  Launcher 

.  B-Z 

None  ' 

B-5 

Elevating  system  —  coil  springs  (total  rate  of  14,000  lb/in.) 
and  little  damping. 

B-6 

Elevating  system  --  coil  springs  and  much  damping. 

B-7 

Elevating  system  --  same  as  B-6,  stiffened  supporting 
structure.  _ 

B-8 

Elevating  system  --  stiff  rods  and  much  damping,  stiffened 
supporting  structure. 

B-9 

Elevating  system  --  ring  springs  (self  damping). 

B- 10 

Elevating  system  --  coil  spring  (low  boost  charge,  720  fps 
muzzle  velocity). 
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All  the  burst  response  curves  are  shown  in  Appendix  J  for  each 
of  the  significant  structural  configurations.  These  response  curves  indicate 
that  the  structure,  although  theoretically  symmetric  about  the  vertical  plane, 
has  an  unsymmetric  response.  The  difference  in  stiffness  between  right  and 
left  elevating  systems  (as  determined  by  measurements)  is  probably  the  major 
cause  of  this  unsymmetrical  response.  Unsymmetrical  supporting  structure 
deflections  were  noted  during  the  measurements  also,  but  these  were  much 
smaller.  These  differences  in  response  from  one  side  of  the  launcher  to  the 
other  increases  the  difficulty  in  analyzing  these  curves  especially  when 
comparing  them  to  theoretical,  planar  response  curves.  Although  these 
unsymmetrical  deflections  can  increase  dUpersion,  we  will  discuss  the 
response  curves  as  though  they  were  identical. 

In  order  to  judge  relative  accuracy  of  these  experimental  response 
curves  let  us  use  a  simple  criterion.  Since  the  recoiling  assembly  is  relatively 
rigid  in  the  plane  of  motion  we  are  considering,  its  motion  relative  to  the 
cradle  due  to  its  own  flexibilities  will  be  at  a  higher  frequency  and  should  damp 
out  more  quickly  than  the  launcher  response.  Therefore,  prior  to  each  sub¬ 
sequent  shot  of  a  burst  the  only  residual  motion  of  significant  magnitude  is  the 
gross  launcher  motion.  The  amount  of  residual  gross  launcher  motion  should 
then  indicate  the  amount  of  dispersion,  since  if  each  round  fires  with  identical 
initial  conditions,  the  response  of  the  launcher  should  be  the  same,  neglecting 
recoil  force  variations  due  to  weight  changes. 

Fig.  J-5  shows  the  cradle  rotation  for  burst  B-5,  fired  with  a 
softened  and  lightly  damped  elevating  system.  This  burst  had  the  greatest 
nonsymmetrical  response  in  comparison  to  all  the  other  firings.  Thi3  could 
be  due  to  additional  nonsymmetry  from  the  dampers.  In  discussing  the 
general  response,  we  will  still  neglect  the  difference  in  the  curves.  The 
cradle  oscillation  for  this  burst  is  a  well  damped  response;  the  maximum  slope 
of  the  response  prior  to  the  second  shot  is  near  0.2  radians/sec.  The  time  of 
shot  ejection  was  not  recorded;  however,  for  maximum  boost  rounds  (muzzle 
velocity  of  1000  fps),  this  time  is  about  20  msec  following  ignition.  The 
ignition  time  is  marked,  and  we  see  that  the  slope  or  velocity  is  near  zero. 

The  trail  strain,  or  supporting  structure  response  in  Fig.  J-6  indicates  ~some 
trail  hop  by  the  relatively  high  frequency  oscillations  about  zero,  immediately 
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following  the  large  initial  response.  This  response  is  also  well  damped. 

Fig,  J-7  shows  the  cradle  rotation  for  burst  B-6  which  had  approximately 
15  times  more  damping  in  the  elevating  system.  As  expected,  the  cradle 
response  decays  more  rapidly,  but  not  as  much  as  was  expected.  The  trail 
strain,  Fig,  J-8,  still  shows  evidence  of  rebound  of  the  structure.  The 
maximum  cradle  velocity  of  the  last  oscillation  prior  to  the  second  shot  was 
reduced  to  0.08  radians/sec.  At  the  e&timated  time  of  shot  ejection  the  cradle 
response  is  very  unsymmetrical  where  one  side  has  a  positive  velocity  near 
0.5 radians /sec  and  the  other  has  a  negative  velocity  near  0.08  radians/sec. 
Theoretically,  the  actual  rotational  velocity  in  the  vertical  plane  is  the 
average,  however,  since  motion  existed  out  of  the  plane,  this  number  may 
be  meaningless.  The  important  fact  is  that  the  overall  response  damped 
out  more  quickly. 

When  the  supporting  structure  was  stiffened  with  the  box  clamps 
and  trunnion  to  trail  struts,  burst  B>7,  the  launcher  response  changed 
significantly.  Fig.  J-9  shows  the  cradle  response  and  Fig.  J-19  the 
corresponding  trail  strain.  The  cradle  rotation  essentially  remained  positive 
at  all  times  as  did  the  trail  strain.  In  addition,  the  magnitude  of  the  cradle 
rotation  was  significantly  lower;  the  peak  trail  strain  was  slightly  lower. 

The  most  important  change,  however,  is  the  shape  of  the  cradle  response  just 
prior  to  the  second  shot;  the  slope  remains  close  to  zero  with  only  small 
variations.  This  indicates  small  cradle  velocity  and  improved  accuracy. 

Another  modification  was  made  to  the  launcher  by  replacing  the 
soft  coil  springs  in  the  elevating  rods  by  relatively  rigid  members  in 
burst  B-8.  The  dampers  and  supporting  structure  stiffeners  were  retained. 
The  response  of  this  system  {Fig.  J- 1 1)  was  very  similar  to  that  of  the 
previous  one  with  possibly  a  small  increase  in  the  residual  cradle  oscillations. 
This  would  seem  to  indicate  that  elevating  rod  flexibility  is  not  a  significant 
parameter.  Another  possibility  exists,  however,  and  is  more  likely  the 
situation.  As  indicated  in  the  section  describing  the  launcher,  significant 
elevating  system  flexibilities  occur  from  sources  other  than  the  rods  them¬ 
selves.  For  this  reason,  when  the  rods  were  stiffened,  the  system  flexibility 
was  not  significantly  changed.  This  is.  due  to  the  fact  that  for  a  series  system 
of  flexibilities,  the  overall  flexibility  is  always  less  than  the  smallest  of  the 
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component  flexibilities.  An  interesting  cradle  response  occurred  during  this 
burst  (three  rounds).  The  cradle  rotation  had  a  larger  mean  residual  value 
between  each  of  the  succeeding  rounds;  in  other  words,  the  response  had  a 
tendency  to  climb.  Although  the  cradle  velocity  was  insignificant  prior  to 
each  shot,  the  displacement  was  increasing.  This  leads  to  an  undesirable 
effect  upon  accuracy. 

Burst  B-9,  Fig.  J-13  and  J-14,  shows  the  result  of  using  ring 
springs  in  the  elevating  system;  ring  springs  have  inherent  damping  since 
their  deflection  is  accompanied  by  sliding  upon  each  other.  These  were 
designed  to  furnish  near  the  maximum  possible  damping  for  the  stiffness 
equivalent  to  the  coil  springs.  Apparently,  this  amount  is  insufficient  to 
furnish  a  good  response. 

The  final  firing  of  a  modified  launcher,  burst  B-lO,  consisted  of 
a  two-round  burst  with  a  zone  7  boost  charge  (muzzle  velocity  of  720  fps). 

The  only  change  in  the  launcher  structure  was  the  addition  of  coil  springs 
into  the  elevating  system;  a  lower  zone  of  boost  ammunition  was  necessary 
to  prevent  bottoming  of  the  springs.  The  cradle  response  (Fig.  J-15)  is 
very  similar  in  shape  to  the  response  of  burst  B-5.  The  newer  response  has 
a  generally  lower  magnitude  because  of  the  lower  projectile  boost.’  It  is  very 
interesting  to  note  that  the  response  is  still  very  well  damped  even  though  no 
damping  was  added.  This  is  pointed  out  again  and  discussed  when  the 
response  of  the  unmodified  launcher  is  discussed. 

A  very  important  characteristic  of  this  launcher  is  illustrated  in 
burst  B-10,  i.  e.  ,  the  variation  in  rate  independent  of  the  response.  The  second 
round  fired  at  approximately  0.7  seconds  after  the  first,  compared  to  less  than 
0.  6  seconds  in  the  other  bursts.  Although  the  variation  in  this  burst  was  later 
discovered  to  be  caused  by  mechanical  interference  in  the  indexing  cams,  other 
changes  can  produce  the  same  effect.  For  example,  changing  the  bypass  valve 
in  the  counterrecoil  flow  can  vary  the  duration  of  counterrecoil  and  therefore 
the  firing  rate.  Also,  variations  in  recuperator  preload  (which  is  sensitive 
-  to  ambient  temperature)  can  also  cause  rate  changes.  The  point  is  that  a 
negative  cradle  Velocity  prior  to  ignition  of  the  rounds  such  as  in  burst  B-5, 
cannot  be  depended  upon  to  decrease  the  value  of  the  initial  positive  response 
following  ignition.  This  is  because  firing  time,  relative  to  the  cradle 
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oscillation,  would  be  very  difficult  to  control.  The  proposed  alternative, 
sought  here,  is  to  reduce  the  cradle  oscillations  to  a  minimum  and  rely  on 
consistency.  This  was  implicit  in  the  simple  criterion  previously  stated. 

One  last  point  to  be  mentioned  about  the  response  of  burst  B-10 
is  the  very  slow  initial  rise  time  of  the  cradle  on  both  shots.  This  is  in 
comparison  to  all  the  other  bursts.  It  is  due  mainly  to  the  slower  rise  time 
of  the  recoil  force  for  the  smaller  boost.  The  type  of  recoil  system  used  has 
this  characteristic  feature. 

We  will  now  discuss  how  these  response  curves  differ  from  those 
of  the  launcher  in  its  as -designed  condition.  Fig.  J-l  shows  the  cradle 
response  during  a  6 -round  burst,  B-2.  In  comparison  to  the  launcher  with  a 
more  flexible  elevating  system,  the  response  is  similar  in  that  the  negative 
cradle  rotation  has  appreciable  magnitude;  the  residual  oscillation,  however, 
has  a  very  small  magnitude  as  well  as  negligible  variation.  (Because  of  the 
perfectly  flat,  residual  response,  the  consideration  of  instrument  error  must 
be  made.  However,  two  other  6-round  bursts  yielded  the  same  result.) 

It  appears  that  the  as -designed  structural  configuration,  by  the  simple  cri¬ 
terion  set  up,  should  produce  the  best  burst  accuracy.  To  a  limited  extent 
this  is  true,  but  there  are  two  important  considerations  which  must  be  men¬ 
tioned.  The  first  consideration  is  the  reason  why  the  as -designed  structure 
responded  the  way  it  did.  If  we  look  at  the  corresponding  trail  strain  curves 
(Eig.  J-2  and  J-3)  we  can  see  that  because  of  the  high  frequency  oscillation 
about  zero,  the  rear  of  the  trails  must  have  lifted  from  the  ground.  This  is 
due  simply  to  the  elastic  rebound,  or  overshoot,  of  the  structure  caused  by 
the  recoil  force  characteristics.  This  means  that  the  supporting  structure 
rotated  forward  about  the  ball  joint  and  compressed  the  two  support  buffers 
between  the  lower  box  and  the  firing  base.  The  amount  of  trail  hop  is  depen¬ 
dent  upon  the  following  quantities: 

1,  Shape  and  magnitude  of  the  recoil  force, 

2.  Drop  off  time  of  the  recoil  force  relative  to  the  natural  period 
of  the  system. 

Variations  in  the  structural  stiffness  or  the  inertia  of  the  launche.  can  alter 
the  occurrence  of  hop.  This  was  shown  when  the  supporting  structure  was 
stiffened  --  the  hop  was  absent. 
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The  second  consideration  is  that  while  the  supporting  structure  is 
tipped  forward,  the  launcher  can  turn  about  a  vertical  axis  through  the  ball 
joint.  This  allows  the  launcher  to  lose  its  aimed  orientation  and  cause 
horizontal  dispersion.  Since  a  method  exists  to  eliminate  this  rotation  during 
hop  (discussed  in  the  next  section),  choice  can  be  made  to  either  eliminate 
the  hop,  by  a  stiffer  carriage  for  example,  or  eliminate  the  rotation  during 
hop.  Either  is  an  acceptable  method  to  insure  good  horizontal  dispersion. 
Other  considerations  may  enter  into  a  choice  of  whether  or  not  trail  hop 
should  be  eliminated,  such  as  the  safety  of  the  crew. 

C .  Horizontal  Dispersion 

Early  tests  of  the  XM70E1  tipping  parts  indicated  that  horizontal 
dispersion  was  larger  than  vertical  and  followed  a  left-right-left  sequence  of 
impact  points  in  burst  firings.  No  analytical  attack  was  made  on  the  disper¬ 
sion  problem  when  it  was  first  recognized  because  at  that  time  no  logical 
theory  of  the  behavior  could  be  constructed.  In  particular,  the  powder 
couple  did  not  offer  a  logical  explanation  because  it  does  not  change  direction 
after  each  shot. 

A  test  program  was  planned  to  investigate  the  sequential 
horizontal  dispersion.  The  weapon  was  prepared  and  the  tests  were  conducted 
with  the  extensive  cooperation  of  the  Rock  Island  Arsenal.  It  was  possible  to 
conduct  the  tests  on  the  force  mount  at  Rock  Island  Arsenal  because  the 
dispersion  evidenced  itself  with  a  high  degree  of  repeatability  during  the 
proof  tests  of  both  Prototype  No.  2  and  Prototype  No.  3. 

The  first  group  of  tests  showed  that  any  round  fired  from  a  breech 
tube  in  the  left  cluster  would  fall  5  in.  or  more  to  the  left  of  the  vertical  bore- 
sight  line  and  that  any  round  fired  from  the  right  would  fall  5  in.  or  more  to 
the  right,  at  an  85-yd  range,  for  all  conditions  tried.  The  boresight  spot  on 
the  target  coincided  before  and  after  each  firing. 

Instrumentation  changes  were  made,  and  the  weapon  was  prepared 
for  a  second  series  of  tests.  The  tests  showed  no  substantial  improvement 
with  slugs  more  accurately  centered  in  the  breech  tubes  or  with  disconnected 
igniter  tubes,  but  some  improvement  was  obtained  with  the  clusters  shimmed 
tight  between  the  front  and  back  plates.  The  instrumentation  indicated  that 
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the  hub  of  the  cluster  being  fired  was  not  being  properly  loaded  with  the 
inertia  force  of  the  cluster. 

Instrumentation  changes  were  nra  de,  and  the  weapon  was  prepared 
for  a  third  series  of  tests.  These  tests  showed  that  a  round  could  be  fired 
accurately  from  either  cluster  if  that  cluster  were  drawi  back  tight  against 
the  back  plate  when  the  round  was  in  the  firing  position.  This  was  accomplished 
by  preloading  the  cluster  thrust  bearing  located  in  the  back  plate  of  the  recoil¬ 
ing  assembly  framework.  Indexing  during  firing  was  then  impossible  and  was 
prevented  by  removing  the  indexing  cam  followers. 

At  this  time,  it  was  realized  that  the  breech  tube  being  fired  was 
being  driven  back  against  the  back  plate  and  was  carrying  the  connected  cluster 
rearward  with  it.  This  resulted  in  the  application  of  a  large  moment  to  the 
rear  plate  since  the  inertia  force  of  the  non-firing  cluster  was  being  applied 
at  its  thrust  bearing. 

A  final  series  of  tests  was  conducted  with  the  breech  tubes 
permitted  a  few  thousandths  of  an  inch  axial  movement  in  the  retaining 
spiders.  This  allowed  the  tube  being  fired  to  move  rearward  without  carrying 
its  cluster  with  it.  The  horizontal  dispersion  was  very  substantially  reduced. 
Several  6-rd  bursts  at  zone  10  (1000  fps  muzzle  velocity)  were  fired  through 
targets  at  85-yd  range.  The  target  pattern  was  nearly  rectangular  with 
overall  dimensions  of  approximately  5  in.  vertical  by  10  in.  horizontal. 

The  remaining  horizontal  dispersion  is  still  sequential. 

Calculations  show  that  because  the  breech  tube  being  fired  is  driven  rearward 
primarily  by  the  powder -gas  pressure  acting  on  a  belled  portion  at  the  front 
of  the  tube,  its  mass  is  lost  from  the  cluster  assembly  and  the  loss  of  that 
mass  then  produces  the  sequential  horizontal  dispersion.  Figure  6  and  7  show 
schematically  how  the  cluster  inertia  forces  are  distributed  throughout  the 
recoiling  framework. 

A  series  of  firing  tests— ^  was  conducted  with  Prototype  No.  3 
during  September  and  October,  I960  at  Redstone  Arsenal  to  determine  the 
effect  of  increased  rigidity  between  the  front  and  rear  plates  (i.  e.  rigidity 
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of  the  recoiling  assembly  framework)  on  the  sequential  horizontal  dispersion. 

A  calculation  showed  that  sequential  horizontal  dispersion  would 
be  very  small  if  the  recoiling  assembly  framework  was  very  rigid.  The 
problem  then  was  to  find  a  lightweight  structure  that  would  greatly  increase 
the  stiffness  of  the  framework.  It  was  decided  that  the  most  efficient 
structure  would  be  a  thin  metal  skin  or  shell  enclosing  the  entire  framework. 

The  firing  program  was  set  up  to  fire  6 -round  bursts  with  and 
without  the  shell.  To  measure  shell  effectiveness,  two  quantities  were 
recorded:  (1)  the  bending  strain  of  the  firing  tube  near  its  base  in  the  hori¬ 
zontal  plane,  and  (2)  the  target  patterns.  Six  bursts  were  fired,  three  with 
and  three  without  the  shell.  (The  launcher  also  contained  additional  stiffening 
fixtures  used  in  the  previous  Fort  Sheridan  firings,  consisting  of  carriage 
clamps  at  the  front  and  rear  corners  of  the  box  sections,  trunnion -to -trail 
struts.) 

The  records  of  bending  strain  in  the  firing  tube  without  the  shell 
exhibited  the  alternating  direction  of  initial  bending  moment.  With  the  shell, 
the  initial  direction  of  the  bending  moment  was  consistently  in  one  direction, 
which  can  he  logically  explained  by  the  rifling  torque;  see  Fig.  8.  The  target 
patterns,  however,  did  not  show  any  reduction  in  horizontal  dispersion.  It 
was  decided  that  the  rigid  body  rotation  of  the  launcher,  due  to  the  trail  pads 
shifting  along  the  ground,  was  hindering  the  effectiveness  of  the  shell. 

A  portion  of  the  test  was  thus  repeated  in  an  attempt  to  obtain 
less  horizontal  dispersion.  This  time  the  trails  were  restrained  from 
rotating  about  the  ball  joint  in  the  plane  of  the  ground  by  placing  a  pinned 
connection  between  a  stake  and  a  trail  on  each  side  of  the  launcher.  This 
permitted  small  rearward  and  upward  trail  motions.  In  this  condition,  four 
6-round  bursts  were  fired,  two  with  and  two  without  the  shell.  The  resulting 
target  patterns  are  shown  in  Fig,  9.  The  two  bursts  with  the  shell  are  rounds 
No.  93  -  98  and  99  -  104;  the  two  bursts  without  the  shell  are  No.  106  -  111 
and  112  -  117.  The  firings  with  the  shell  had  a  total  horizontal  dispersion  of 
approximately  1. 5  and  1. 9  mils,  while  the  two  without  the  shell  had  a  total 
horizontaldisper  sion  of  3.  7  and  2.  9  mils.  Because  of  the  single -hole  targets, 
it  is  impossible  to  calculate  the  standard  deviation,  but  it  can  be  inferred 
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that,  for  the  two  bursts  with  the  shell,  the  standard  deviation  is  less  than 
1  mil. 

The  results  of  these  tests  indicate  that  the  shell  offers  a  reduc¬ 
tion  of  horizontal  dispersion  of  nearly  100%.  (The  trial  ties  showed  a  reduc¬ 
tion  of  200  to  300%.)  The  reliability  or  consistency  is  fairly  satisfactory  since 
these  targets  represent  four  consecutive  bursts  fired  from  Prototype  No.  3 
under  identical  conditions.  But  in  order  to  obtain  these  results,  it  must  be 
emphasized  that  the  trails  were  restrained  from  rotating  about  a  vertical 
axis  through  the  ball  joint. 
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IV.  THEORETICAL  ANALYSES 

The  theoretical  analyses  of  the  launcher  dynamics  consist  of 
two  specific  mathematical  models  to  determine  launcher  motion,  and  three 
supplemental  studies  which  used  very  simple  dynamic  models  in  order  to 
directly  study  the  effects  of  certain  parameter  variations.  The  two  models 
which  were  developed  to  calculate  the  launcher  motion  are: 

1.  A  planar,  nonlinear,  3-degree -of-freedom  system, 

2.  A  three-dimensional,  nonlinear,  n-degree -of-freedom  system. 

The  three  supplemental  studies  include: 

1.  A  linear,  2-degree-of-freedom  analysis  to  study  the  effectiveness  of 
various  locations  for  auxiliary  damper  locations, 

2.  A  third-order,  linear  single-degree-of-freedom  model  to  determine 
optimum  elevating  system  structural  parameters. 

3.  A  linear,  second-order  system  to  determine  the  parameters  which 
bring  about  an  optimum  response  decay. 

The  derivation,  assumptions,  and  some  results  of  these  models 
and  studies  are  presented  in  the  appendices;  a  discussion  of  their  development, 
applications  and  results  are  presented  in  this  section. 

The  derivation  and  discussion  of  a  mathematical  accuracy  criterion 
is  also  included  in  this  section.  This  criterion  is  based  upon  the  sensitivity  of 
the  rocket  to  the  different  launcher  motions  studied  in  the  mathematical  models 


Planar,  Nonlinear,  3-Degree-of-Freedom  Model 


Prior  to  the  initiation  of  this  Long-Range  Study  Program,  some 
dynamic  studies  were  conducted  which  aided  the  formulation  of  this  model. 
One  study,  in  particular,  was  a  similar,  planar,  3-degree -of-freedom  model 
with  the  following  degrees  of  freedom: 

1.  Rigid  body  rotation  of  the  launcher  about  the  trail  ends;  base  hop, 

2.  Recoil  motion, 

3.  Elastically  restrained  rotation  about  the  ball  joint. 


The  last  type  of  motion  was  assumed  to  the  first  mode  response  of  the 
relatively  flexible  supporting  structure,  with  a  rigid  tipping  mass  and  rigid 
elevating  system. 
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Subsequent  experience,  both  from  the  use  of  the  above  model  and 
from  actual  launcher  firings,  provided  valuable  information  for  making 
assumptions  for  the  new  model-  For  example,  instrumented  firings  and  the 
output  of  the  above  model  both  showed  that  rotation  about  the  trail  ends  was 
very  small.  It  was  therefore  neglected  in  the  new  model.  Recoil  motion 
calculated  in  the  above  model  agreed  fairly  well  with  expe  riment  except  during 
counterrecoil.  The  difference  was  traced  to  seal  friction  and  indexing  forces. 

The  former  results  in  longer  counterrecoil  durations,  while  the  latter  brings 
about  two  changes.  Indexing  forcer  cause  small  changes  in  the  recoil  motion 
but  more  important,  they  significantly  affect  the  cradle  response.  Consequently, 
they  were  both  included  in  the  new  model.  Inability  to  obtain  a  good  correlation 
of  the  experimental  response  with  the  output  of  the  prior  model  led  to  the  con¬ 
clusion  that  the  elevating  system  was  flexible.  Again,  this  flexibility  was  there¬ 
fore  included.  The  assumption  of  the  first  mode  response  of  the  supporting 
structure  was  supported  by  previous  experimental  firings,  and  was  reconfirmed 
with  the  dynamic  analysis  described  earlier;  this  assumption  was  retained. 

The  three  degrees  of  freedom  of  the  model  developed  on  this  program  correspond 


to: 

1.  Plastically  restrained  rotation  of  the  supporting  structure  about  the 
ball  joint,  j. 

2.  Recoil  motion  relative  to  the  cradle,  u, 

3.  Rotation  of  the  tipping  parts  about  the  trunnion  relative  to  the 
supporting  structure,  (p  £• 

The  last  motion  is  a  consequence  of  the  flexible  elevating  system.  Although 
no  auxiliary  damping  in  the  elevating  system  was  included  in  the  as -designed 
launcher,  linear  damping  terms  were  placed  in  the  equations  for  dampers  in 
parallel  with  the  elevating  rod3  and  also  a  damper  between  the  lower  box  and 
the  ground. 


The  frictional  forces  were  added  in  the  form  of  a  constant  frictional 
force  due  to  the  preload  plus  a  term  proportional  to  the  recoil  pressure.  This 
corresponds  to  the  type  of  seals  employed.  The  equation  of  motion  of  the 
recoiling  mass  in  the  u  direction  also  includes  terms  representing  the  constant 
force,  constant-stopping  distance  recoil  system  used  in  the  launcher. 
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The  final  equations  of  motion  consisted  of  three,  second-order, 
nonlinear  differential  equations.  In  order  to  solve  these  without  making  any 
restrictive,  simplifying  assumptions,  they  were  programmed  for  ARF's 
Univac  1105  in  the  Univac  Scientific  Exchange  (USE)  language.  This  program 
is  described  in  Appendix  D.  Solution  of  the  three,  second-order  differential 
equations  was  effected  by  using  the  Runge-Kutta  numerical  integration 
technique.  Launcher  parameters  such  as  damping,  stiffness,  recoil  rod 
shape,  initial  boost  charge  force,  etc.  were  all  included  as  input  data;  a 
sample  output  sheet  is  shown  in  Fig.  0-2. 

The  purpose  of  the  computer  solutions  of  launcher  motion  can  be 
grouped  into  three  divisions  by  considering  their  purpose.  The  first  group  of 
solutions  was  aimed  at  matching  the  experimental  launcher  response.  The 
second  was  to  verify  the  existence  of  an  optimum  response  and  determine  the 
optimum  parameters.  (This  is  described  later  when  the  optimum  system  is 
discussed.)  The  final  group  of  computer  runs  was  to  show  the  results  of 
various  changes,  e.g. ,  low  zone  response,  efforts  of  elevation  changes,  and 
delayed  recoil. 

1.  Correlation  of  Theory  with  Experiment 

The  success  of  correlation  between  computed  and  experimental 
response  curves  varied  and  depended  upon  launcher  parameters;  in  general 
the  agreement  was  better  when  less  negative  carriage  rotation  occurred. 

The  difficulty  was  due  to  the  fact  that  the  original  model  did  not  represent 
the  nonlinear  stiffness  and  damping  displayed  by  the  launcher  supporting 
structure.  This  was  subsequently  corrected  and  closer  correlation  obtained, 
though  not  entirely  satisfactory. 

Figures  10  through  21  show  typical  information  obtained  from  the 
computer  program.  Figure  18  shows  a  full  burst  response.  The  numerical 
input  corresponded  to  either  measured  or  calculated  values  of  physical 
parameters  of  Prototype  No.  3  Launcher.  Certain  parameters  were  varied, 
however,  to  obtain  the  correlation.  All  of  the  parameters  not  designated  on 
the  curves  are  listed  on  the  computer  program  in  Appendix  D. 

The  two  curves  shown  superimposed  in  Fig.  23  are  the  cradle 
response  from  Fig.  15  and  the  experimental  cradle  response  from  firing 
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R-38,  listed  in  Table  1-1.  This  experimental  response  was  from  a  single  shot 

with  identical  launcher  conditions  to  burst  B-8,  Fig,  J-ll.  The  supporting 

structure  stiffness  measured  for  this  condition  was  5.  5  x  10  in.  -lb/ radian. 

The  numerical  value  of  carriage  stiffness  used  in  the  computer  program- was 

4.  5  x  107  in.  -lb/ radian  and  is  significantly  lower.  The  fact  that  this  lower 

value  was  necessary  to  match  the  experimental  response  is  probably  due  to 

the  flexibility  added  to  the  system  by  the  soil  in  the  actual  firing.  (The 

structural  stiffness  measurements  were  made  on  concrete.  )  The  measured 

5 

elevating  system  stiffness  for  shot  R-38  was  approximately  4.  5  x  10  compared 

A 

to  1.4  x  10  used  in  the  computer.  The  difference,  though  very  large,  is 
probably  due  to  the  unmeasured  elevating  system  flexibilities  discussed  earlier. 
(Burst  B-7  had  an  elevating  system  stiffness  of  1.4  x  10  lb/ in.  but  had  a  very 
similar  response  to  burst  B-8,  even  though  the  latter  supposedly  had  a  much 
iitiffer  elevating  system.  This  further  supports  the  conclusion  that  the  high, 
measured  stiffness  value  is  not  correct.) 

A  further  comparison  can  be  made  between  the  computed  carriage 
rotation  (Fig.  13)  and  the  trail  strain  from  burst  B-8,  Fig.  J-ll.  These 
curves  agree  in  less  detail,  but  their  main  characteristics  are  the  same.  The 
computed  curve  is  generally  smoother  since  it  does  not  contain  any  higher 
mode  response  characteristics. 

Figures  26  and  27  show  experimentally  obtained  strain  from  the 
elevating  screw  and  damper  supports.  These  curves  are  tensile  and  com¬ 
pressive  strain  in  circular,  pinned  rods  and  are  directly  proportional  to  the 
forces.  They  compare  reasonably  well  with  the  spring  force  and  damper 
force  shown  in  Fig.  18  and  19  respectively. 

Because  of  the  similarity  of  the  response  of  burst  B-7  and  B-8, 

we  see  that  the  computed  response  also  agrees  reasonably  well  with  B-7.  In 

4 

this  case  the  actual  launcher  had  an  elevating  stiffness  of  1.4  x  10  .  Burst 
B-6  ,  Fig.  J-7,  reveals  a  large  negative  dip  after  the  initial  positive  displace¬ 
ment,  which  did  not  occur  in  burst  B-7  or  B-8.  The  only  difference  is  that 
B-6  did  not  have  a  stiffened  supporting  structure.  Reduction  of  only  the 
supporting  structure  stiffness  in  the  computer  model  produced  the  response 
:.n  Fig.  24.  The  expected  correspondence  does  not  occur;  the  cradle  response 
remains  essentially  positive.  The  negative  supporting  structure  response, 
however,  doubled  compared  to  the  previous  theoretical  response  which  indi- 
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Fig.  10  TYPICAL  COMPUTER  SOLUTION  --  POWDER  GAS  FORCE 


11  TYPICAL  COMPUTER  SOLUTION  --  RECOIL  FORCE 
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Fig.  12  TYPICAL  COMPUTER  SOLUTION  --  RECOIL  DISPLACEMENT 
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ig.  14  TYPICAL  COMPUTER  SOLUTION  --  CARRIAGE  ROTATIONAL.  VELOCITY 
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Fig.  15  TYPICAL  COMPUTER  SOLUTION  --  CRADLE  ROTATION 


Fig,  16  TYPICAL  COMPUTER  SOLUTION  --  CRADLE -CARRIAGE 
RELATIVE  ROTATION 


TYPICAL  COMPUTER  SOLUTION  --  ELEVATING  SYSTEM  SPRING  FORCE 


19  TYPICAL  COMPUTER  SOLUTION  --  ELEVATING  SYSTEM  DAMPER  FORCE 


TYPICAL  COMPUTER  SOLUTION  --  AXIAL  CAM  FORCE 
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Fig.  23  COMPARISON  OF  COMPUTED  AND  EXPERIMENTAL  RESPONSE  CURVES 


cates  that  the  nonlinear  negative  supporting  structure  characteristics  are 
very  significant  here.  To  investigate  the  effects  of  a  bilinear  stiffness  the  *• 
computer  program  was  modified.  Figure  25  shows  the  cradle  response 
using  a  bilinear  supporting  structure  stiffness.  Although  the  agreement  is 
much  better,  it  could  probably  be  improved  with  slightly  different  values  of 
stiffness  and  also  the  additional  simulation  of  a  bilinear  damping  coefficient. 
The  effect  of  these  nonlinearities  has  been  demonstrated,  however,  even 
though  perfect  correlation  was  not  obtained. 


2.  Effects  of  Elevation  Changes 


It  was  important  to  investigate  the  effects  of  elevation  changes 
upon  the  response  of  the  launcher.  This  was  done  easily,  because  the 
elevation  was  an  input  to  the  computer  program.  Figures  28  and  29  show 
the  cradle  response  of  the  launcher  with  structural  parameters  identical 
to  those  in  the  previously  computed  curves.  These  responses  are  at  20°  and 
60°  elevation,  respectively,  flhe  previous  ones  were  at  (45°).  There  are  some 
definite  changes  in  the  details  of  the  response,  especially  at  the  time  imme¬ 


diately  following  the  ignition  of  the  rounds.  At  20  elevation,  the  initial 


rise  time  (velocity  at  shot  ejection)  is  significantly  higher  than  at  45°,  and 


at  45°  the  initial  response  is  higher  than  that  corresponding  to  60°.  There 


is  a  very  simple  reason  for  this.  The  geometry  change  brings  about  a  change 
in  the  inertial  moment  of  the  tipping  assembly.  To  demonstrate  this,  Fig.  30 
illustrates  that  as  the  supporting  structure  deflects  due  to  the  application  of 
the  recoil  force,  the  trunnion  moves  in  a  direction  perpendicular  to  a  line 
through  the  trunnion  and  ball  joint.  At  20°  elevation,  this  motion  causes  the 
tipping  assembly  to  react  with  an  initial  counterclockwise  rotation  which  adds 
to  the  normal  response.  At  60°,  however,  the  initial  tendency  of  the  tipping 
assembly  is  to  rotate  clockwise  which  opposes  the  normal  counterclockwise 


rotation.  The  result  here  is  a  delay  in  the  cradle  response,  which  occurs 

o 

at  every  round  of  the  burst.  The  response  at  45  elevation  is  similar  to 


60  ,  but  the  delay  is  smaller  in  magnitude.  We  now  have  a  situation  where 

increases  in  elevation  bring  about  greater  accuracy  --  an  inherent  charac¬ 
teristic  of  this  launcher.  This  characteristic  has  less  importance  when 
delayed  recoil  is  used  and  is  apparent  when  delayed  recoil  is  discussed  in 
Section  IV -E. 
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EXPERIMENTAL  ELEVATING  SYSTEM  DAMPER  SUPPORT  STRAIN 


Fig.  29  CRADLE  RESPONSE  AT  60°  ELEVATION 

AHMOUP  PESEAUCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


3.  Low  Boost  Firings 

This  launcher  has  the  capability  of  firing  rockets  with  various 
muzzle  velocities  (zones  of  ammunition)  simply  by  varying  the  amount  of 
gunpowder.  The  powder-gas  force  curve  corresponding  to  a  muzzle  velocity 
of  380  fps  was  used  in  the  computer  program.  Figure  31  shows  the  cradle 
response  for  60°  elevation  at  this  low  boost.  (The  high  value  of  elevation 
was  necessary  to  insure  a  large  enough  recoil  distance  to  cause  automatic 
indexing.)  As  would  be  expected,  the  rise  time  of  the  initial  cradle  response 
is  much  smaller  than  for  the  high  boost  rounds  and  was  due  to  the  slower 
rise  time  in  the  recoil  force.  The  amplitude  of  the  residual  motion  prior  to 
successive  shots  is  also  less,  but  it  is  still  significant.  The  variation  in 
the  residual  motion,  however,  is  very  similar  to  the  high  zone  bursts  since 
the  counterrecoil  forces  developed  are.  a  function  of  the  recuperator  force 
and  the  indexing  forces.  These  are  both  relatively  independent  of  the  boost, 
because  of  the  constant  stopping  distance  recoil  system.  The  maximum 
recoil  distance  for  this  burst  was  near  20  inches  compared  to  24  at  higher 
zones.  As  a  result  of  the  slower  rise  times  of  the  recoil  force,  we  see  that 
it  is  much  more  possible  to  have  a  negative  cradle  velocity  at  shot  ejection. 

In  addition,  delayed  recoil  will  have  less  of  an  effect  in  changing  the  response. 

4.  Importance  of  Coriolis  Acceleration 

Because  of  the  high  recoil  velocity  (400  in.  /sec)  combined  with 
the  rotational  velocity  of  the  tipping  assembly,  a  question  arose  concerning 
the  importance  of  the  Coriolis  acceleration.  This  acceleration  is  proportional 
to  the  product  of  these  two  velocities.  For  this  mathematical  model,  because 
of  computer  solution,  there  was  no  particular  reason  to  neglect  this  term;  in 
other  situations,  this  might  lead  to  a  convenient  simplification  in  the  equa¬ 
tions  of  motion.  For  this  reason,  Coriolis' effect  upon  the  cradle  response 

was  investigated.  Because  the  term  is  proportional  to  the  product  of  the 
■■  ■  .  •  • 

velocities  u  and  (  ^)*  a  set  of  parameters  was  selected  where  these 

velocities  were  high.  The  recoil  velocity  u  is  almost  independent  of  the 
structural  parameters,  so  the  selection  was  reduced  to  finding  parameters 
which  increased  (  tjp  ^  .  This  is  done  by  increasing  the  stiffness  of 

the  elevating  system.  Figures  32  and  33  show  the  cradle  response  with  and 
without  the  Coriolis  terms.  The  overall  response  shape  has  not  changed 
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while  the  peak  response  has  a  difference  of  7%.  This  indicates  that  the 
Coriolis  acceleration  is  of  second  order  importance  with  respect  to  the  shape  of  the 
cradle  response,  though  the  change  in  magnitude  may  be  significant.  This 
serves  to  indicate  that  some  of  the  simplified  models  used  in  the  supplementary 
analyses  may  be  a  reasonable  representation  of  the  launcher  for  predicting 
large  scale  changes  in  the  response. 

5.  Delayed  Rocoil 

The  recoil  system  on  the  XM70  Launcher  makes  use  of  a  variable 
area  orifice  which  was  designed  to  stop  the  recoiling  assembly  in  a  constant 
distance  with  a  constant  force.  The  orifice  area  varies  as  a  function  of  the 
recoil  distance;  this  was  effected  by  pulling  variable  diame  ter  rods  through 
a  fixed  area  orifice  plate.  A  design  suggestion  was  made  to  modify  the  recoil 
system  in  such  a  way  that  the  recoil  force  is  zero  (or  of  negligible  magnitude) 
until  the  rocket  has  left  the  firing  tube.  In  this  way  the  carriage  deflections 
will  be  zero  (for  the  first  shot  of  a  burst),  the  velocity  will  be  small,  and  the 
launcher  accuracy  will  be  near  perfect.  Certain  problems  arise  with  delayed 
recoil;  the  magnitude  of  the  recoil  force  increases  due  to  the  shorter  distance 
over  which  the  recoil  force  acts,  because  of  the  constant  stopping  distance 
system.  A  question  also  arose  about  the  effect  of  the  delay  on  shots  following 
the  first. 

The  3-dcgree-of-freedom  mathematical  model  was  used  to  study 
the  effect  of  delayed  recoil  since  the  recoil  rod  shape  is  entered  into  the  model 
as  a  tabular  function  of  the  recoil  displacement.  Inserting  a  rod  shape  into 
this  mathematical  model  corresponds  identically  to  inserting  new  recoil  rods 
into  an  actual  launcher. 

The  first  rod  shape  tried  was  a  standard  rod  with  a  decreased  dia¬ 
meter  (increase  in  orifice  area)  for  the  first  11-1/2  in.  of  recoil  motion.  The 
first  rocket  leaves  the  launcher  at  about  3  to  4  in.  of  recoil;  subsequent  shots 
are  fired  by  an  adjustable  cam,  9  in.  out  of  battery.  Because  of  the  residual 
counterrecoil  velocity  during  all  shot  ignitions  except  the  first,  these  rockets 
leave  the  launcher  at  a  recoil  distance  near  11-1/2  in.  Figure  34  shows  the 
cradle  rotation  and  the  recoil  force  much  larger  than  obtained  for  the  standard 
rod.  Figure  5  shows  a  peak  force  near  20, 000  lb  compared  with  the  new  force 
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Fig.  34  CRADLE  RESPONSE  FOR  DELAYED  RECOIL  ROD  CUTBACK  11 


near  30,000  lb,  The  later  shots,  however,  are  near  the  20,000  lb  peak 
reflecting  the  fact  that  of  the  11-1/2  in,  cutback,  only  about  2-1/2  in.  of  it 
affect  the  later  shots,  compared  to  the  entire  11-1/2  for  the  first  shot. 
Comparison  of  the  response  to  that  shown  in  Fig.  20  (optimum  elevating  system 
parameters)  shows  that  delayed  recoil  cut  the  cradle  rotation  by  about  50%  at 
shot  ejection  (approximately  20  msec  after  t  =  0)  and  lowered  the  velocity 
by  approximately  30%.  The  displacement  effect  for  later  shots  is  almost 
negligible  because  it  is  fairly  large  to  begin  with.  Close  inspection  of  the 
curves  shows  that  the  velocity  reduction  is  significant  for  the  second  shot. 

The  third  rocket,  however,  left  the  launcher  when  the  cradle  had  a  negative 
velocity  (for  the  regular  control  rod).  The  use  of  the  cutback  rods  changed 
the  burst  characteristics  such  that  the  third  shot  fired  prior  to  a  region  of 
positive  cradle  velocity.  Wo  then  realize  the  situation  that  delayed  recoil  does 
not  aid  the  accuracy  if  the  cradle  would  normally  have  a  negative  velocity  at 
the  time  a  round  is  ignited.  The  normal  launcher  response  tends  to  lessen 
this  negative  velocity  with  a  desirable  result.  This  leads  us  to  the  conclusion 
that  in  order  to  gain  the  benefits  of  delayed  recoil,  the  residual  cradle  motion 
should  be  as  small  as  possible.  In  this  situation,  delayed  recoil  is  definitely 
useful.  V 

In  order  to  decrease  the  high,  first-round  recoil  force  a  rod  was 
tried  with  two  diameter  cutbacks,  the  first  from  0  to  3  in.  and  the  second 
from  8  to  11-1/2  in.  It  was  expected  that  the  extra  raised  portion  between 
3  and  8  in.  would  help  decrease  the  first  force  level.  Figure  35  shows  the 
result;  the  peak  force  is  reduced  from  30,000  to  26,  000  lb.  The  peak  cradle 
rotation  was  also  reduced.  The  later  rounds  were  unaffected  because  the 
section  of  the  rod  over  which  they  operate  was  unchanged. 

We  saw  earlier  that  firings  at  60°  elevation  have  less  initial  cradle 
rotation  due  to  the  geometry  of  the  elevating  system  and  tipping  assembly  mass 
center.  A  computed  response  for  a  60°  firing  using  a  cutback  recoil  rods  is 
shown  in  Fig.  36.  In  this  situation,  the  cradle  velocity  is  nowhere  near  that 
of  the  cradle  for  the  ordinary  launcher.  In  the  same  respect,  the  response 
for  20°  elevation  (with  a  cutback  rod)  is  probably  near  the  ordinary  response. 
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Pig.  35  CRADLE  RESPONSE  FOR  DELAYED  RECOIL,  NOTCHED  ROD 


Pig.  36  CRADLE  RESPONSE  FOR  DELAYED  RECOIL  ROD  CUTBACK  11 


B.  Supplementary  Analyses 
1..  Third-Order  System 

A  third-order  linear  system  was  developed  as  a  mathematical 
model  and  used  to  find  the  stiffness  and  damping  which  should  be  incorporated 
into  the  elevating  system  of  the  launcher.  This  model  was  obtained  by 
simplification  of  the  launcher  representation  shown  in  Fig.  C-l.  The 
simplifications  were: 

1.  The  inertia  of  the  carriage  and  trails,  IQi  was  assumed  negligible 
compared  to  the  moment  of  inertia  of  the  tipping  assembly  --  a 
ratio  of  approximately  115:8000. 

2.  Nonlinear  coupling  terms  between  the  motion  of  tipping  assembly 
and  supporting  structure  were  neglected. 

3.  The  recoiling  motion  was  neglected. 

4.  The  recoil  force  was  treated  as  the  primary  disturbance  and  was 
approximated  by  an  impulse. 

These  assumptions  are  considered  valid  because  the  resulting  solution  is 
not  intended  to  represent  the  actual  launcher  motion,  but  rather  to  indicate 
the  variation  of  the  response  behavior  for  a  wide  range  of  parameter 
variations.  The  linear  analog  of  the  resulting  model  is  sl\pwn  in  Fig.  E-l. 

The  responses  of  this  third-order  system  were  determined  by 
using  the  response  plots  contained  it:  "Dynamic  Response  Plots  and  Design 
Charts  for  Third-Order  Linear  Systems"  —  .  This  report  contains  the 
solution  of  such  systems  in  the  form  of  dimensionless  response  plots  for  forty- 
eight  different  combinations  of  parameters.  These  plots  were  obtained  from 
evaluation  of  Equation  E-5,  as  a  function  of  time,  on  a  digital  computer. 

The  design  charts  in  the  report  1/  show  areas  of  three  types  of 

responses . 

1.  Parameter  combinations  which  produce  a  normal  damped  response. 

2.  Parameter  combinations  which  produce  an  unstable  response. 

3.  Parameters  which  result  in  three  real  roots  to  Eq.  E-3. 

In  order  to  establish  the  undesirability  of  the  third  type  of  response,  a  set  of 

1/  Meyfarth,  P,  ,  Dynamic  Response  Plots  and  Design  Charts  for  Third- 

Order  Linear  Systems.  Research  ivlemo  No.  &.  M.  7401-3,  Massachusetts 
Institute  of  Technology,  1958. 
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parameters  which  produce  this  wan  used  in  the  computer  solution  for  the 
3-degree-of-freedom  system.  The  resulting  cradle  response  was  continuous 
ly  increasing  (non-oscillating),  as  shown  in  Fig.  E-3,  for  150  lb-sec/in. 
damping  and  100  lb/in.  stiffness. 

Despite  the  various  assumptions,  the  response  curves  obtained  were 
a  fair  representation  of  the  actual  launcher  motion  for  a  single  shot.  (Com¬ 
pare  Figures  E-2  and  E-4.)  The  optimum  response  for  this  system  was  de¬ 
fined  as  one  with  both  a  low  first  peak  displcaement  and  a  low  ratio  of  second 
peak  to  first  peak  displacement.  Three  response  plots  which  satisfy  this 
definition  are  shown  in  Figure  E-2.  The  corresponding  optimum  elevating 
system  parameters  are  a  stiffness  of  1.0  x  10  to  1.4  x  10  lb/in.  and  damp¬ 
ing  of  360  to  540  lb-sec/in. 

2.  An  Approximate  Dynamic  Analysis  of  the  Minimum  Weapon  Weight 

Required  for  a  Specified  Accuracy 

A  relatively  simple  mathematical  model  of  the  XM70  and  a  simple 
dynamic  performance  requirement  for  accuracy  was  formulated  in  order  to 
evaluate  several  basic  design  parameters  including  the  minimum  total  weight 
of  the  weapon.  Mathematically,  the  procedure  was  one  of  finding  a  second 
order  system  which  had  a  specified  response  decay.  The  details  of  the 
analysis  are  contained  in  Appendix  B  of  this  report. 

For  the  purposes  of  this  analysis,  the  complete  weapon  was  represent¬ 
ed  by  a  linear,  forced,  damped,  single -degree-of-freedom oscillator.  The 
displacement  of  the  oscillator  represented  angular  rotation  of  firing  tube 
in  a  vertical  plane.  This  representation  was  shown  to  be  useful— ■  ■ when  the 
stiffness  is  taken  to  be  the  total  stiffness  of  the  weapon  structure  and  the  mass 
is  a  computed  dynamic  equivalent.  Values  of  the  damping  coefficient  were 
chosen  to  be  consistent  with  experimental  data  of  the  weapon's  response  to 
firing  loads.  The  force  applied  to  the  oscillator  to  simulate  firing  was  a 
rectangular  pulse  whose  amplitude  and  duration  were  chosen  to  represent 
the  recoil  force  of  the  actual  weapon. 

As  a  criterion  for  adequate  dynamic  performance,  it  was  assumed 
y  See  ARF  Project  K130,  Bimonthly  Report  No.  11,  Appendix  A. 

AS  MODS  UESEARCH  POUNDATION  OP  ILLINOIS  INSTITUTE  OP  TECHNOLOGY 

-67-  . 


that  the  amplitude  of  firing  tube  oscillations  (in  a  vertical  plane)  must  be  less 
than  some  small  value,  which  is  considered  to  produce  negligible  error  at 
the  target,  prior  to  firing  of  each  round  of  a  burst. 

Under  the  assumptions  of  the  analysis,  a  closed  form  solution  can  be 
obtained  to  the  equation  of  motion  of  the  weapon.  The  solutions  consists  of 
the  displacement  of  the  firing  tube  as  a  function  of  time  and  the  stiffness, 
damping  and  mass  of  the  weapon  but  what  is  ultimately  desired  from  the 
analysis  is  the  time  required  for  adquate  damping  of  the  response  as  a  func¬ 
tion  of  the  weapon  weight. 

The  weight  of  the  weapon  was  explicity  introdxccd  into  the  solution  by 
writing  expressions  for  the  stiffness,  damping  coefficient,  and  equivalent 
mass  of  the  weapon  as  functions  of  the  weight  of  the  major  subassemblies 
of  the  XM70.  These  expressions  were  obtained  by  assuming  linear  rela¬ 
tionships  and  using  the  XM70  Prototype  Model  No.  1  weights  for  reference. 
Using  the  solution  containing  these  expressions,  it  was  possible  to  obtain  a 
time -displacement  curve,  and  hence  the  time  required  for  adequate  damping, 
for  a  particular  total  weapon  weight.  Several  solutions  must  be  obtained  in 
order  to  obtain  a  plot  of  required  damping  time  versus  total  weight  from 
which  the  minimum  weight  can  be  obtained.  The  analysis  procedure  is 
facilitated  by  using  the  envelope  of  the  response  curves  rather  than  the 
curves  themselves. 

This  analysis  formed  the  basis  for  the  computation  of  the  optimum 
recoiling  assembly  weight  for  the  XM70  which  is  reported  upon  in  ARF 
Project  Kl 30  Bimonthly  Report  No.  8  Appendix  B  .  In  addition,  this 
analysis  was  used  to  estimate  the  minimum  weight  of  a  weapon  having  a 
small  recoiling  assembly  weight  (500  lb)  and  a  delayed  recoil  force.  The 
minimum  weight  of  such  a  weapon  was  found  to  be  only  1340  lb.,  but  the 
resulting  ground  reaction  was  too  large  and  the  calculated  carriage  weight 
was  inadequate  from  a  strength  standpoint.  These  calculations  led  to  the 
development  of  an  auxiliary  strength  criterion  for  the  carriage  and  hence 
to  the  procedure  referred  to  above  for  calculating  optimum  recoiling 
assembly  weight. 

3.  Linear,  2-Degree-of-Freedom  Model 

A  classical,  2-degree-of-freedom  system  was  another  system  studied 
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to  obtain  information  about  the  best  location  for  auxiliary  damping  of  the 
cradle  motion.  Because  only  qualitative  information  was  desired,  the  lack 
of  representation  of  system  nonlinearities  is  not  considered  important. 

Fig.  A-l  shows  the  system  which  was  studied  where  m ^  represented  the 
tipping  assembly,  andm^,  the  supporting  structure.  The  three  damper 
locations  which  were  considered  correspond  to  the  following: 

1.  c. ;  damping  between  the  supporting  structure  and  the  ground  (the 
ground  was  assumed  to  be  capable  of  furnishing  a  tensile  load). 

2.  c  ;  damping  between  the  supporting  structure  and  tipping  assembly. 

u 

3.  Cj?  damping  between  the  tipping  assembly  and  the  ground. 

The  equations  of  motion  were  solved  by  digital  computation  where  the  input 
disturbance,  Fj,  was  a  tabular  representation  of  an  actual  recoil  force. 

All  three  damper  locations  are  physically  attainable,  though  the  second  is 
certainly  the  most  convenient  location.  Equivalent  damping  between  the 
three  locations  was  determined  by  estimating  the  amount  of  forcible  damping 

for  the  same  weight  contribution  to  the  launcher  from  the  respective  dampers. 
* 

As  expected,  damping  directly  between  the  cradle  and  ground  was 
most  effective,  of  the  three  locations.  However,  damping  between  the  tipping 
assembly  and  supporting  structure  was  relatively  effective,  especially  com¬ 
pared  to  damping  between  the  supporting  structure  and  ground.  Combinations 
of  the  latter  two,  Cj  and  c^,,  was  better,  but  on  the  estimated  weight  basis, 
the  improvement  was  not  significant.  Since  c^-damping  is  the  most  diffi¬ 
cult  to  attain,  the  most  promising  damper  location  was  concluded  to  be 
relative  to  the  supporting  structure  and  tipping  assembly. 

C.  Three-Dimensional,  n-Degree-of-Freedom  Model 

This  mathematical  model  (called  Mathematical  Model  VI  in  progress 
reports)  was  developed  with  a  view  toward  the  incorporation  of  several  of  the 
modes  of  motion  which  were  not  included  in  earlier  models .  It  contains  all 
of  the  freedom  of  the  previous  models  and,  in  addition,  permits  the  investi¬ 
gation  of  the  following  types  of  behavior. 

1.  Ground  Flexibility 

ResuTts'of  sxirrre  of  the  -experiments  showed  that  an  elastic  and  plastic 
deformation  of  the  earth  underneath  the  ball  joint  and  the  trail^patha  can 
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sizably  affect  the  burst  accuracy  of  the  weapon.  The  new  model,  conse¬ 
quently,  will  allow  the  introduction  of  an  arbitrary  force -displacement 
relationship  at  the  support  points. 

2.  Sidewise  Rotation  at  the  Ball  Joint 

Incorporation  of  this  freedom  necessitated  that  a  full  three-dimon- 
sional  representation  be  utilized.  This  representation  should  permit 
simulation  of  a  portion  of  the  lateral  dispersive  action  associated  with  the 
weapon.  ■  • 

3.  Left-Right  Sequencing 

This  is  also  a  three-dimensional  effect;  it  arises  because  the 
inertial  unbalance  in  the  recoiling  assembly  cannot  be  eliminated.  Left- 
right  sequencing  is  discussed  in  detail  in  another  portion  of  this  report. 

Its  incorporation  into  this  model  was  somewhat  simplified  through  the  in¬ 
terpretation  of  the  recoiling  assembly  framework  flexibility  in  the  form  of 
a  single  torsional  spring. 

4.  Flexible  Barrel 

The  barrel  flexibility  was  Considered  to  be  very  important  to  the 
overall  burst  dispersion  of  this  weapon.  In  particular,  it  would  appear  that 
phasing  between  the  barrel  vibration  and  the  time  of  firing  of  the  various 
rounds  is  quite  significant.  This  barrel  vibration  is  even  more  important 
because  it  is  evidently  coupled  with  the  rotation  of  the  weapon  about  the  ball 
joint  to  Coriolis  force  terms  and  gyroscopic  moments  which  arise  due  to  the 
moving  projectile. 

5 .  Moving  Projectile 

The  flexible  barrel  is  excited  by  the  traveling  load  induced  through 
a  projectile  as  well  as  through  the  gross  forces  which  arise  through  firing. 
This  operation  motivated  the  inclusion  in  the  model  of  a  representation  for 
the  coupled  motion  of  a  traveling  mass  and  the  barrel  proper.  This 
representation  is  so  organized  that  an  experimentally  determined  time- 
displacement  relationship  can  be  introduced  to  describe  the  axial  projectile 
behavior.  It  is  thus  possible  to  determine  the  instantaneous  muzzle  position 
orientation  and  velocity  at  the  time  of  projectile  exit.  It  is  felt  that  this 
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information  should  bo  of  aid  in  the  development  of  more  complete  accuracy 
criteria  for  rocket  launchers  of  this  type. 

This  3-dimenaional  mathematical  model  was  not  completed  until 
near  the  close  of  the  present  project.  The  equations  were  programmed  for 
the  UNIVAC  1105  digital  computer  at  Armour  Research  Foundation,  but  the 
program  was  not  run  on  the  machine  because  it  was  felt  that  the  time  and 
efforts  which  remained  on  the  project  could  be  more  efficiently  applied  to 
other  facets  of  launcher  behavior.  The  equations  for  this  model  are  pre¬ 
sented  in  Appendix  F  of  the  present  report.  These  equations  are  quite 
complex;  their  full  solution  should  prove  to  be  time  consuming.  Conversely, 
specific  items  of  interest  to  future  rocket  launcher  programs  can  certainly 
be  obtained  through  specilization  of  these  very  general  equations. 

D.  Accuracy  Criterion 

The  main  direction  of  effort  on  this  program  has  been  to  study  the 
response  of  the  XM70  Launcher  structure  in  the  vertical  plane  due  to 
burst  firing  loads.  The  intention  was  to  find  that  combination  of  launcher 
parameters,  e.  g.  .  weight,  stiffness,  moment  of  inertia,  mass  distribution, 
recoil  force,  etc.  ,  which  would  optimize  or  at  least  improve  the  response. 

In  order  to  evaluate  the  effects  of  parameter  changes,  a  simple  accuracy 
criterion  was  used.  This  criterion  was:  the  most  accurate  burst  has  the 
smallest  residual  angular  cradle  displacement  and  velocity  prior  to  each 
succeeding  round.  This  was  not  expressed  mathematically,  but  was  applied 
simply  by  examining  graphs  of  both  computer  and  experimental  responses. 

A  mathematical  accuracy  criterion  is  presented  here  that  permits  a  more 
accurate  basi3  for  judgment.  This  criterion  is  a  relationship  between  the 
launcher  response  and  the  points  of  projectile  impact.  General  application 
of  this  criterion  is  limited  only  by  the  fact  that  it  employs  the  motion  charac¬ 
teristics  of  a  short-burning -time,  ballistic  rocket.  Because  the  development 
of  this  criterion  came  late  in  this  program,  its  application  was  not  completed 
It  is  presented  here,  however,  because  of  its  generality. 

The  significant  sources  of  rocket  dispersion  (of  a  rocket  with  a 
relatively  short  burning  time)  are: 

1.  Initial  cross  spin,  i.  e.  ,  angular  yaw  (pitch)  velocity  when  the 
rocket  leaves  the  launcher, 
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2.  Launch  mal-point,  i.e.,  initial  yaw  (pitch)  when  the  rocket 
leaves  the  launcher, 

3.  Dynamic  unbalance  of  rocket, 

4.  Cross  winds  during  burning, 

5.  Total  rocket  impulse  variations, 

6.  Rocket  spin  at  burnout, 

7.  Thrust  malalignment, 

Only  the  first  three  of  these  are  influenced  by  the  launcher.  The 
particular  launcher  motion  that  can  contribute  significantly  to  dispersion  by 
influencing  these  three  quantities  is  the  cradle  motion  in  combination  with 
the  firing  tube  motion  relative  to  the  cradle.  The  cradle  angular  displacement 
and  angular  velocity  in  the  vertical  plane  are  considered  here.  At  present, 
we  are  limited  to  consider  only  these  because  all  completed  mathematical 
models  assume  a  rigid  firing  tube  and  plane  motion  only.  Because  of  these 
two  limitations,  dynamic  unbalance  must  also  be  neglected.  In  addition, 
during  the  transition  period  when  the  rocket  is  partially  out  of  the  firing 
tube,  the  rocket  is  assumed  to  have  the  orientation  of  the  tip  of  the  tube- 

The  purpose  of  this  criterion  is  to  permit  a  comparison  of  launcher 
responses  to  burst  firings  (or  a  group  of  single  shot  firings  with  the  same 
launcher)  and  a  subsequent  selection  of  the  best  response;  with  this  in  mind, 
the  following  characteristics  can  be  imposed  upon  the  criterion: 

1.  Only  relative  changes  in  dispersion  or  accuracy  between  burst 
responses  need  be  indicated  by  the  criterion. 

2.  In  the  case  of  computed  responses,  an  indication  of  accuracy 
improvement  by  the  mathematical  output  of  the  criterion  must 
correspond  to  an  actual  accuracy  improvement.  (This  can  be 
checked  experimentally  or  with  a  complete,  ballistics  mathe¬ 
matical  model).  This  correspondence  will  not  be  required  to 
be  linear,  however. 

3.  The  criterion  will  yield  a  numerical  answer  for  each  burst  and 
will  be  zero  for  a  perfect  burst.  (A  perfect  burst  is  attained 
when  all  rounds  impact  at  the  same  point. ) 

4.  The  criterion  must  allow  for  the  different  contributions  to 
dispersion  from  unit  deviations  in  initial  yaw  and  initial  cross 
.spin. 

5.  Each  round  of  a  burst  is  to  be  considered  separately,  and  its 
sequential  position  will  be  neglected. 
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These  are  the  criterion  characteristics  that  are  inferred  prior  to 
the  development  of  the  criterion  simply  by  considering  its  purpose  and 
application.  The  implications  of  each  of  these  are:  Item  1  states  that  the 
criterion  output  for  a  single  burst  will  not  indicate  anything  about  the 
accuracy  of  that  burst,  when  compared  to  the  criterion  output  of  another 
burst;  however,  it  will  indicate  which  of  the  two  is  better,  i.  e. ,  more 
accurate,  Item  2  states  that  the  indication  of  which  of  two  bursts  was  more 
accurate  must  correspond  to  the  indication  of  relative  accuracy  of  an  actual 
launcher  firing  of  two  bursts  under  the  same  conditions;  the  criterion  must 
agree  with  reality.  In  order  that  the  criterion  contain  the  characteristics 
implied  by  Item  3,  the  deviations  of  both  displacement  and  velocity  of  the 
cradle  will  be  measured  with  respect  to  their  mean.  Item  4  states  that 
the  criterion  must  account  for  the  fact  that  a  unit  deviation  in  angular 
cradle  displacement  does  not  necessarily  produce  the  same  dispersion  as 
a  unit  deviation  in  angular  cradle  velocity;  this  is  accomplished  by  using 
what  are  known  as  "unit  effects",  coefficients  determined  by  the  magnitude 
of  dispersion  due  to  a  unit  deviation  in  each  quantity.  The  last  item.  No.  5 
merely  states  that  no  priority  is  given  to  any  of  the  six  rounds  of  a  burst. 

The  combination  of  the  above  characteristics  narrows  the  choice 
of  mathematical  expressions  describing  this  accuracy  criterion  to  a  point 
where  an  expression  can  be  deduced: 

L  -  Uv  £  1a,  I  *  V 


k  =  1 
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U 
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U 


the  numerical  output  from  the  above  expression 
indicating  relative  burst  accuracy, 

unit  effect  for  mal  -point  or  initial  yaw, 

unit  effect  for  initial  cross  spin, 

deviation  from  the  mean  in  initial  yaw, 

deviation  from  the  mean  in  initial  cross  spin. 

index. 


The  two  deviations  are  written  as: 
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0k  -  0, 


where  (p  and  (pm  are  the  initial  yaw  and  cross  spin,  respectively. 

In  this  criterion,  the  smaller  the  value  of  L,  the  better  the  burst.  Other 
criteria  may  be  deduced  that  conform  to  the  given  conditions,  but  for  small 
deviations  they  will  reduce  to  the  above  expression. 


This  criterion  was  to  be  incorporated  directly  into  the  computer  pro¬ 
gram  of  the  3-degree-of-freedom  mathematical  model.  A  6-rd  burst  com¬ 
putation  would  then  automatically  give  the  quantity.  L.  The  unit  effects 
are  obtained  from  solutions  of  the  equations  of  motion  of  the  rocket;  the 
unit  effects  for  the  XM54  rocket  were  tabulated  in  previous  progress  reports 
of  this  project.  ~ 


—  For  more  information  of  the  unit  effects  of  rocket  dispersion,  see: 
Bi-Monthly  Report  No.  18,  ARF  Project  D124,  January  22,  I960. 
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V.  CONCLUSIONS 

All  the  result*  derived  from  this  study  are  based  upon  considera¬ 
tions  of  cradle  motion  near  the  time  of  .shot  ejection  of  each  round  of  a 
burst.  There  was  no  intention  of  neglecting  the  effects  of  firing  tube 
flexibilities,  but  this  effort  represents  a  first  step  to  increase  the  weapon 
accuracy.  There  was  no  reason  to  suspect  that  flexible  firing  tube  motions 
may  be  simply  additive  to  cradle  motion  because  of  the  coupling  which 
exists.  However,  the  correlation  obtained  between  theory  and  experiment, 
despite  neglecting  this  coupling  theoretically,  indicates  that  the  effects 
upon  gross  cradle  rotation  are  small.  This  points  to  the  conclusion  that 
the  accuracy  criterion  based  upon  cradle  motion  was  valid,  i.  e. ,  a  reduc¬ 
tion  in  the  residual  cradle  motion  prior  to  the  successive  rounds  of  a  burst 
will  improve  accuracy. 

The  residual  cradle  oscillations  of  the  as-designed  launcher  (Proto¬ 
type  No.  3)  were  seen  experimentally  to  be  small.  The  only  disturbing 
factor  was  that  the  as  designed  structure  was  subjected  to  elastic  rebond, 
or  trail  hop;  this  permitted  serious,  gross  horizontal  dispersion  to  occur. 
This  horizontal  dispersion  was  eliminated  by  simple  ties  restraining  the 
trails  from  moving  in  a  horizontal  plane.  An  alternative  was  shown  to 
exist  where  trail  hop  was  eliminated  --  this  was  shown  to  be  accomplished 
by  stiffening  the  supporting  structure.  With  this  stiffened  structure,  the 
ties  at  the  trails  are  unnecessary. 

When  trail  hop  does  not  occur,  an  optimum  combination  of  elevating 
system  structural  parameters  was  found.  This  was  based  upon  a  study  which 
showed  that  the  best  location  (based  upon  effectiveness  per  weight)  at  which 
auxiliary  damping  should  be  added  to  the  launcher  is  between  the  tipping 
assembly  and  the  supporting  structure.  The  optimum  elevating  system 
had  the  following  structural  parameters: 

4 

(total  stiffness)  =  1.4x10  lb/in. 

c_,  (Damping  coefficient)  =  400  lb-sec/in. 

This  combination  caused  a  minimum  of  residual  cradle  oscillation  between 
shots  of  a  burst,  and  is  independent  of  the  stiffness  of  the  supporting 
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structure.  The  possibility  of  an  unstable  response  during  a  burst  was 
demonstrated  for  very  low  elevating  system  stiffness.  It  should  be 
mentioned  that  due  to  the  damping  of  the  hydraulic  recoil  system  and 
possibly  soil  interaction,  the  launcher,  as  designed,  displayed  much 
inherent  damping. 

Elevation  changes  were  shown  theoretically  to  influence  the  accuracy 
because  of  the  corresponding  geometry  changes.  The  overall  offset  was  to 
increase  accuracy  with  elevation.  Changes  in  accuracy  were  also  observed 
due  to  difference  in  ammunition  boost  charge.  The  difference,  howover, 
was  much  less  than  might  be  expected  because  of  the  large  reduction  in  the 
recoil  force  at  the  lower  zones.  This  was  because  the  forces  developed 
during  counterrecoil,  from  the  recuperator  and  indexing  system,  are  of 
the  same  order  of  magnitude  as  during  the  high  boost  bursts. 

Delayed  recoil  is  effective  in  reducing  the  velocity  and  displacement 
of  the  cradle  response  for  single  shots.  For  bursts,  the  effectiveness  is 
questionable  since  it  will  improve  accuracy  when  the  ignition  of  the  round 
occurs  simultaneously  with  positive  cradle  velocity,  but  it  does  not  help 
when  ignition  occurs  during  a  duration  of  a  negative  cradle  velocity.  Rods 
cut  back  to  delay  the  recoil  force  for  the  first  shot  only  will  probably  yield 
the  best  overall  accuracy  improvement.  This  improves  single  shot  accuracy 
while  it  eliminates  the  possibility  of  detracting  during  later  shots. 

When  considering  horizontal  dispersion,  the  shell  which  stiffens 
the  recoiling  assembly  framework  reduces  the  dispersion  by  more  than 
100%.  In  addition,  when  trail  hop  is  present,  horizontal  ties  between  the 
trails  and  the  ground  are  necessary  to  achieve  good  accuracy.  They  may 
also  improve  dispersion  in  the  absence  of  hop  if  a  tendency  exists  for  the 
trail  pads  to  slide. 

Finally,  inclusion  of  Coriolis  acceleration  was  shown  to  be  significant 
in  the  calculation  of  peak  response  magnitude,  but  neglecting  it  did  not  alter 
the  overall  response  shape. 

In  order  to  completely  investigate  launcher  contributions  to  disper¬ 
sion,  two  additional  efforts  are  recommended.  These  are  a  completion  of 
the  development  of  the  three-dimensional  mathematical  model  begun  in  this 
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project,  and  a  comprehensive  experimental  program  utilizing  the  stiffen* 
irg  fixtures  developed.  This  experimental  program  should  include  firing 
tube  motion  instrumentation  as  well  as  rocket  instrumentation.  It  is  felt 
that  this  combined  effort  can  produce  a  weapon  whose  accuracy  in  burst 
fire  exceeds  that  of  present  single  shot  artillery. 
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APPENDIX  A 

A  STUDY  OF  THE  RELATIVE  EFFECTIVENESS 
OF  VARIOUS  DAMPER  LOCATIONS 

N.  P.  Pearson,  R.  M.  Brach,  and  R.  H.  Van  Beek 

I.  INTRODUCTION 

This  appendix  concerns  the  theoretical  possibility  of  damping 
the  motions  of  the  supporting  structure  and  tipping  parts  of  an  XM70-type 
launcher  between  consecutive  shots  of  a  burst.  Ideal  viscous  dampers  are 
considered  to  be  located  between  the  ground  and  supporting  structure,  and/or 
between  the  supporting  structure  and  the  tipping  parts,  and/or  between  the 
ground  and  the  tipping  parts. 

II.  THE  MATHEMATICAL  MODEL 

Figure  A-l  is  a  schematic  drawing  of  an  XM70-type  launcher 
simplified  and  idealized  to  the  point  where  it  has  only  two  degrees -of  freedom 
The  general  procedure  for  forming  such  a  discrete  system  from  an  elastic 
system  consists  of  lumping  equivalent  masses  at  points  where  the  stiffness 
is  known  and/or  loads  are  applied  and  determining  the  values  of  the  equiv¬ 
alent  masses  in  such  a  way  that  the  actual  and  reduced  systems  have  the 
same  kinetic  and  potential  energy  during  their  first-mode  oscillation.—^ 

The  accuracy  of  the  technique  depends  upon  the  nature  of  the  forcing  func¬ 
tion  and  upon  the  location  where  the  motion  is  desired.  Deflection  magni¬ 
tudes  at  any  given  instant  are  not  expected  to  be  accurate  because  impor¬ 
tant  features,  such  as  the  indexing  forces,  are  not  considered.  However,  the 
relative  effectiveness  of  the  three  damper  locations  is  expected  to  be  correct¬ 
ly  represented.  Clearly,  in  this  particular  case,  if  the  calculations  showed 
the  motions  to  be  thoroughly  damped  after  some  relatively  long  time,  say  tQ, 
and  if  the  actual  weapon  were  equipped  with  analogous  dampers,  the  motion 
of  the  actual  carriage  or  cradle  would  be  thoroughly  damped  after  also. 

1.  Jacobsen,  L,  S,  ,  Ayre,  R.  S.  ,  Engineering  Vibrations,  McGraw  Hill 
Book  Co.  ,  Inc.  ,  New  York,  1958. 
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The  equations  of  motion  of  the  system  shown  in  Fig,  A-l  are; 

-  ■- ‘  2  ~  <•  l»  '  *1  A-‘ 


and 

(Ie  +  m  J})  \  =  k2  {02  -  0L)  +/F(t)  -  kj  0j  -  Cl  0j  A-2 

or,  letting  +  c£  {  02  -  ^  ) 

'e  *  mT^  "  ‘l 

and 

i  F(t)>  M(t) 

then 

IjiJj  -  k2(02  -  0^+ kj  0!  +  ct  +  c2  {  02  -  ^  M  (t)  A-3 


+  k2(02  -01>+  c2(  02  -  0!  )  +  3  02  A -4 

Equations  A-3  and  A-4  have  been  solved  with  a  forcing  function, 

M  (t),  derived  from  an  experimental  record  of  recoil  pressure.  Figure  A-2 

shows  the  forcing  function  used  for  all  the  solutions  given  in  this  appendix. 

The  value  of  the  function  is  just  the  recoil  pressure  multiplied  by  the 

appropriate  area  and  moment  arm.  Figures  A-3  through  A-8  show  the 

calculated  response  for  various  choices  of  the  parameters.  The  stiffness 

values,  and  k^,,  represent  the  values  either  measured  or  calculated 

for  the  existing  weapons.  In  addition  to  the  values  of  damping  and  stiffness, 

each  of  the  graphs  is  labeled  with  the  maximum  value  of  the  quantity 
•  • 

c2  (0i  -  0^)  which  is  the  maximum  force  developed  in  the  damper  between 
the  supporting  structure  and  the  tipping  parts. 
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III.  DISCUSSION  AND  CONCLUSIONS 

Figure  A-3  describe*  the  undamped  response  of  the  system, 
which  is  not  representative  of  the  actual  weapon  response,  primarily  be¬ 
cause  the  actual  weapon  contains  inherent  damping,  The  figure  does  show, 
however,  that  significantly  larger  amplitudes  may  be  obtained  on  shots 
after  the' first.  This  phenomenon  was  observed  experimentally  for  Prototype 
No.  1.  It  is  worthwhile  noting  that,  in  all  of  these  response  plots,  the 
negative  values  of  0^  may  be  taken  as  a  rough  indication  of  the  system's 
tendency  to  produce  trail  hop. 

The  parameters  used  in  the  calculation  represented  by  Fig.  A-4, 

A-5,  and  A-6  are  identical  except  for  the  values  of  damping,  with  c,,  c,, 

5  2  *  * 

and  Cj  alternately  taking  on  the  value  1.25x10  lb-sec  /in.-,  the  other 
two  being  zero  in  each  case.  Note  from  Fig.  A-6  that  damping  between 
the  ground  and  the  tipping  parts,  c^-damping,  gives  the  best  results.  The 
relative  effectiveness  of  C2~damping,  as  shown  by  Fig.  A-5,  and  the 
relative  ineffectiveness  of  c^ -damping,  as  shown  by  Fig.  A-4,  constitute 
the  most  important  result  of  this  study. 

Damping  between  the  tipping  parts  and  the  ground,  c^-damping, 
is  the  most  difficult  type  to  obtain  practically.  Figure  A-7  shows  an  attempt 
to  obtain  the  effectiveness  of  c^-damping  using  only  c^-  and  c^-damping. 

The  results  are  seen  to  be  not  as  good  as  c^-damping  alone.  This  is  partic¬ 
ularly  discouraging  when  the  weight  added  by  the  two  dampers  is  considered. 
Figure  A-8  shows  the  results  of  attempts  to  compensate  for  the  weight  of  the 
two  dampers  by  reducing  the  supporting  structure  stiffness  and,  hence,  its 
weight.  This  performance  is  seen  to  be  slightly  less  satisfactory,  but 
further  conclusions  are  unwarranted  because  of  the  simplicity  of  the  model 
which  does  not  include  indexing  forces  and  elevating  system  geometry. 
Because  these  factors  affect  the  optimum  stiffness -damping  relationship, 
this  relationship  will  be  studied  by  means  of  a  more  elaborate  mathematical 
model. 
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Fig.  A-2  MOMENT  PRODUCED  BY  THE  RECOIL  FORCE 


CALCULATED  RESPONSE  OF  AN  XM70-TYPE  LAUNCHER  WITH  c, -DAMPING  ONLY 


CULATED  RESPONSE  OF  AN  XM70-TYPE  LAUNCHER  WITH  c-,  -  DATING  ONLY 


Fi g  A -7  CALCULATED  RESPONSE  OF  AN  XM70-TVPE  LAUNCHER 


APPENDIX  B 


A  DYNAMIC  ANALYSIS  OF  SOME  BASIC  DESIGN  PARAMETERS 

By  R.  H,  Van  Beek 


I.  INTRODUCTION 

This  appendix  contain*  a  dynamic  analysis  which  was  intended  to 
aid  in  determinimg  some  of  the  basic  design  parameters,  such  as  weight 
and  stiffnesses,  of  weapons  of  the  XM70  type.  This  analysis  uses  a  relative* 
ly  simple  mathematical  model  of  a  weapon  together  with  semiempirical  data 
from  the  XM70  Prototype  Model  No.  1.  A  nomenclature  list  is  at  the  end  of 
this  Appendix. 

II.  ANALYSIS 

It  has  been  shown  that  the  weapon  can  be  represented  by  a  sin- 
gle-degree-of-freedom  dynamic  model  with  the  recoil  force  (rod  pull)  as  a 
forcing  function.  The  assumed  model  and  forcing  function  are  shown  in 

figures  B-l  and  B-2. 

F 

i  F,  lb 


k,  U) 


time,  sec 


Fig.  B-l  MATHEMATICAL 


Fig.  B-2  FORCING 


The  differential  equation  for  the  system  of  figure  B-l  is 


M  =  x  +  cx  +  kx  =  F 


Eq.  (B-l) 


Let  yZ  t0  • 


t  — y  XM 

»  U  —  *  >  ZZ  '  \ 


Eq.  (B-2) 
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The  solution  to  differential  equation  (B-l)  in  terms  of  these 


quantities  is 

xm 


■i^[l  -  e  -yr  (  ainYl/^ 


com  Vf  f  1 


*t)  for  0  *  Of  *  1 


Eq.  (B-3) 

~r'] 


Xl?)  =  e'*/4'*'  [l^sin  yj  J\  2  +a  L2cos  V),/l  - 

for  *  1  Eq.  (B-4) 


Eq.  (B-5) 


and  3 


sin 


Eq.  (B-6) 


For  a  single  degree -of-freedom  torsional  system,  the  same 
solution  holds  with  the  dimensionless  variables  given  by 


x  = 


TT7  •  V  ‘o'  /*"  2 fa  '  r'  S, 

O  O  v/ 


Eq.  (B-7) 


A  typical  displacement  time  plot  is  shown  in  figure  B  -3. 

Adequate  dynamic  performance,  i.e.  acceptable  accuracy,  results 
when  t^  is  less  than  the  time  to  a  subsequent  firing  and  9^  is  a  displacement 
which  produces  negligible  error  at  the  target.  Several  primary  design 
parameters  can  be  logically  determined  by  finding  the  minimum  weight 
weapon  which  will  give  the  required  performance. 
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For  givoiyz  and  ^  the  graph  of  5C  versus  can  be  drawn.  The 
last  intersection  of  the  curve  with  0  =  +_  0C  determines  Repeating  this 
procedure  for  a  number  of  combinations  o t Jul  and  Y[  gives  data  for  a  plot  of 
7^  versus ^  for  various  values  oi^C,  For  a  given^,  a  plot  ofig  versus 
JLL  can  then  be  drawn.  All  points  on  this  curve  are  such  that  the  correspond* 
ingVJ  and^tt  produce  a  displacement  of  ©c  at^=^and  0  <  0c  for'?'>'2£.  All 
points  give  values  otyu  and  >7  for  which  ?£is  less  than  the  value  of  ^  for 
which  the  curve  is  drawn.  Figure  B-4  shows  a  typical  graph  of'fversus^  » 

As  an  approximation,  the  envelope  of  equation  (B-4)  can  be  used. 

It  is  given  by 


,2  .  t  2 

L1  +  L2 


Eq.  (B-8) 


Then 

(0e) c  =  +  LZ2  Eq.  (B-10) 

This  eliminates  the  necessity  of  drawing  the  first  two  of  the  above 
mentioned  graphs,  since,  given  7c  and  9c,  for  ea.chy<  the  corresponding  vf 
can  be  found  numerically  by  trial  and  error.  In  all  cases  where  the  graphs 
of^tt  versus  were  found  by  both  methods  the  difference  between  the  graphs 
was  negligible. 

Let  the  total  weight  of  the  weapon  be  given  by 

Wt  =  WL  +  W2  +  W3  f  WNS  Eq.  (B-ll) 

Assume  that  the  weight  of  the  carriage  is  related  to  its  rotational 
stiffness  by  the  linear  relation 

W  =  <Xk  Eq.  ( B -12) 
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B-4 


Eq.  (B-18) 


III.  RESULTS 

From  this  relationship,  curves  of  ^  versus^  for  various  values 
of  W^,  can  be  plotted.  If  a  curve  of  ^  versus/X  for  the  desired  value  of  1? 
is  plotted  on  the  same  graph,  the  lowest  possible  weight  giving  acceptable 
accuracy  can  be  determined.  The  situation  is  illustrated  in  Fig.  5.  The  * 
constant  ^curve  can  be  thought  of  as  representing  acceptable  performance 
and  the  other  curves  arc  the  constant  weight  curves.  Then  the  lowest  weight 
curve  which  just  intersects  the  performance  curve  gives  the  minimum  total 
weight  and  the  point  of  intersection  determines  the  desired  damping  ratio 
and  .  Then  by  the  use  of  equations  (B-ll),  (B-12),  and  (B-13),  the  car¬ 
riage  stiffness  and  weight  and  the  damper  weight  can  be  calculated.  Fig.  5 
also  shows  the  results  of  applying  this  analysis  with  the  following  input  data: 

Projectile  weight*  50  lb. 

Muzzle  velocity  1000  ft/s cc. 

Firing  rate  20rnd6/min. 

Maximum  Amplitude  of 

Vibration  at  triggering  (0  )  0.0005  rad. 

Maximum  Recoil  Displacement  24  in. 

(Free  Recoil  until  shot  ejection  is  assumed  and 
requires  12  in.  of  recoil) 

Recoiling  Assembly  Weight  500  1b. 


The  inertia-weight  relationships  and  the  carriage  stiffness  -  weight  relation¬ 
ship,  (oC),  of  the  XM  70  were  assumed.  It  will  be  noted  that  the  hypothetical 
weapon  has  a  very  low  firing  rate  and  small  recoiling  assembly  weight  and 
hence  has  a  small  mimimum  weight.  From  the  graph,  the  minimum  total 
weight  which  gives  the  desired  accuracy  is  1340  lb.  Included  in  this  is  a 
total  accessory  or  non-structUral  weight  which  is  assumed  to  be  750  lb. 

This  is  the  same  as  the  XM70  and  includes  everything  except  the  recoiling 
assembly  and  the  carriage  structure.  Using  the  values  of yW  and  Yj  from  the 
graph,  the  damper  weight  is  found  to  be  67.5  lb.  and  the  carriage  weight 

these  quantities  determine  the  forcing  function 
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22.5  1b.  The  carriage  stiffness  required  is  1. 12  x  lO^lb.  in. /rad. 

This  example  design,  although  it  theoretically  satisfies  the  accur¬ 
acy  criterion,  is  impossible  because  a  carriage  weighing  22.5  lb.  could  not 
carry  the  rod  pull.  This  result  led  to  the  development  of  strenghth  criterion 
and  a  method  of  optimizing  recoiling  assembly  weight.  This  latter  method 
was  used  for  the  XM70  Prototype  Model  No.  2  and  is  presented  in  ARF  Pro¬ 
ject  K130  Bi-monthly  Report  No.  Appe  dix  B. 
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IV.  NOMENCLATURE 

WT 

total  weight  of  gun 

WNS 

weight  of  gun  hardware 

W1 

weight  of  recoiling  assembly 

W2 

weight  of  carriage 

W3 

weight  of  damper 

k  torsional  stiffness  of  carriage 

K  spring  stiffness  for  linear  vibration  system 

c  damping  stiffness  for  linear  vibration  system 

M  mass  stiffness  for  linear  vibration  system 

b  =  m 

“  A 

I  effective  moment  of  inertia  of  torsional  vibration  system 


=  /-j—  ,  natural  frequency  of  undamped  system 


ratio  of  damping  to  critical  damping 


7^  =  CO  i  non-dimensional  frequency 


displacement  for  linear  vibration  system 


9  displacement  for  torsional  vibration  system 


C  t  ime 
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B-  10 


Mx 


*0 

Toto2 


non-dimensional  displacement 


length  of  pulse 


amplitude  of  pulse-linear  vibration  system 


amplitude  of  pulse-torsional  vibration  system 


t/t  non-dimonsional  time 
'  o 


moment  of  inertia  of  recoiling  assembly  about  the  base 


moment  of  inertia  of  carriage  about  the  base 
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APPENDIX  C 

EQUATIONS  OF  MOTION  OF  THE  LAUNCHER 
WITH  THREE  DEGREES  OF  FREEDOM 

R.  M.  Brach 


I.  INTRODUCTION 

Thu  following  aquations  apply  to  launchers  with  a  structural 
geometry  corresponding  to  that  of  the  XM70  automatic  rocket  launcher.  It 
is  a  three -degree -of-freedom  analysis  of  gross  launcher  motion  in  the 
vertical  plane.  These  equations  are  to  be  used  to  investigate  the  stability  * 
of  the  system  and  the  influence  on  stability  of  changes  in:  (1)  elevating 
system  stiffness  and  damping,  (2)  carriage  and  trail  stiffness,  (3)  mass 
distribution,  (4)  recoil  force,  and  (5)  indexing  cam  length.  They  can  also 
furnish  an  indication  of  the  effects  due  to  variations  in  the  muzzle  momentum 
of  the  projectile.  However,  because  the  analysis  considers  a  rigid  firing 
tube,  it  cannot  predict  the  precise  projectile  motion  at  the  time  of  muzzle 
ejection. 

The  three  degrees  of  freedom  (see  Fig.  C-l)  are  the  following: 

(1)  the  linear  displacement  of  the  recoiling  mass  relative  to  the  cradle,  u, 

(2)  the  angular  motion  of  the  tipping  parts  relative  to  the  carriage,  0 
and  (3)  the  absolute  angular  motion  of  the  supporting  structure,  (f)  j.  The 
supporting  structure  in  this  analysis  is  taken  to  be  the  entire  structure  that 
supports  the  tipping  parts;  i.  e.  ,  the  supporting  structure  includes  a  rigid, 
ira  ssless,  trunnion  side,  an  inertia,  and  a  massless,  flexible  trail.  The 
inertia  is  an  equivalent  inertia  found  by  letting  the  above  support  contain 
kinetic  and  potential  energy  equal  to  that  of  the  actual  launcher  vibrating  in 
its  static  deflection  shape.  The  supporting  structure  is  assumed  to  rotate 
about  the  ball  joint,  point  0  in  Fig.  C-l;  a  provision  for  damping  the 
rotation  is  included.  The  elevating  system  is  assumed  to  be  a  massless 
spring  and  damper.  The  cradle  and  recoiling  mass  are  both  assumed  to  be 
rigid  bodies. 

The  equations  of  motion  of  the  recoiling  mass  include  those 
describing  the  constant-force,  constant-stopping -distance  recoil  system, 

A9MOUS  BtSEAUCH  FOUNDATION  OP  (lll.NOIS  INSTITUTE  OP  TECHNOLOGY 


. 


used  in  the  XM70,  and  also  the  axial  reactions  of  the  indexing  cam  since 
these  have  been  found  to  be  significant  from  past  experience. 

II.  ANALYSIS 

A.  Kinematics 

If  the  motion  of  the  recoiling  mass  is  described  by  xf  and  yf 
with  respect  to  the  ground  and  by  u  with  respect  to  the  cradle,  we  see  from 
Fig.  C-2  that 

Xf  =  JL  COS  </?♦  0,1  +  u  cos(OC  +  0  l  +  0 2> 

and  C-i 

yr  =  /sin(^}+01)  +usin(CX+0j+02) 

In  these  equations,  JL  is  the  distance  from  the  point  0  to  the  trunnion, 
point  T  ,  fj  in  the  fixed  obtuse  angle  made  with  ground  and  the  trunnion 
sides,  a  i»  the  angle  of  elevation,  0  ^  is  the  angular  motion  of  the  support¬ 
ing  structure  and  0^  is  the  angular  motion  of  the  tipping  parts  with  respect 
to  the  supporting  structure.  The  velocities,  xr  and  y^  ,  and  the  accelera¬ 
tions,  xr  and  y^,  are  found  by  differentiation  to  be 

*r  *  .in  1/^4  0,)  4  u  co«  (  0C+ 0, +02) 

-«(01,02>  •  10(0(40, 402> 
yr=  i0  J  cos  1/1+0,) +  S  sin  (  CX  +  0j  +  0  2) 

+  u  (  01  +  02^  cos  ^  +  0  i  +  02^  C-2 

and 

xr  =  -j?  0  l  sin  (0?  +  0L)  -  J  0J2  cos  ( y3  +  0  t) 

+  u  cos  (  CX  +  0l  +02)  -  U  (  0  J  +  0,)  sin  (  0C+  0  1  +  02^ 

-  u  (  0  j  +  02)  sin  ((CX  +  0j  +  (p  2)  -  u  (  0^  +  02)  sin 

(CX+  0,  +02)  -  U  (  0 1  4  02)2  co,  (C<  +'01+02)  „  ■ 
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yr  3  l  cf)l  COB  (  JJ+pJ  -  JL(piZ  *rn(/3+  +  U  sin 

( a + 0 !  +  0  2)  +  u  ( $  j  +  cj)2)  cos  <a  +  0!  +  0  2> 
cos  (CX+  ^  +  ^2^  +  **  (  ■+  (p2) 

co«(a+0l+02)  -u(('/)1+^)2)2»in{a  +  0l+02)  C"4 

The  vector  eum  of  the  component  and  yf  may  be  resolved  into 
components  in  the  u-  and  v-  directions.  These  components  are  denoted 
as  and  £ ,  respectively.  From  Fig.  C-3,  these  are: 

7}  3  xf  cos  (Jl  +  yr  sin  (J  ^  C-5 

£3  Vr  COi  CJj  *  xf  »in  Uj  I  C-6 

where 

Ui  3  CX+01  +02  .  C-7 

Substituting  Eq.  C-3  and  4  into  Eq.  C-5  and  6  gives  the  following  expression. 
7|  =  -  X  (j)  j  sin  L/2  -  Ji.  (j)^  cos  +  u  -  u  (  Cf)l  +  (p%)2  C-8 

£  3  i$x  cos  U2  -l(plZBinUz*Zn^l  +<^2)  +  u((J1+(^2). 

C-9 

where 

V2  *f$-<X-cpz.  c-10 


B.  Equations  of  Motion 

.  Equations  C-8  and  C-9  are  the  components  of  the  absolute 

acceleration  of  the  recoiling  mass.  The  free  body  diagram  of  the  recoiling 
mass  is  shown  in  Fig.  C-4.  Summing  forces  acting  through  the  mass  center 
we  obtain 

mr  jj  =  -  F  (t)  +'<J)(t)  +  ey*A  j  N| -  m^g  sin  U^.  C-ll 

£  =  N  -  m^g  cos  LJ ^  .  C-12 

In  these  equations  the  following  nomenclature  is  used: 

F  (t)  applied  force 
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T*  (t)  applied  moment  (powder  couple  equal  to  £  F  (t)  where  6"  is  the 
moment  arm  from  the  bore  center  to  the  recoiling  matt  eg. 

(^)  (t)  reactions  of  the  recoil  system  and  cam 

N  normal  reaction  force  between  the  cradle  and  recoiling  mass 

T  moment  reaction  between  the  cradle  and  recoiling  mass 

n  coefficient  of  Coulomb  friction  between  the  cradle  and 

recoiling  mass 

e  e  =  1,  u<0;  e  a  -1,  u»0;  e  =  0,  u  =  0 

The  XM70  launchers  have  very  low  friction  ball  bearings  between  the  cradle 
and  recoiling  parts;  the  force  fj\ n|  is  therefore  negligible  compared  to  the 
seal  friction.  Because  the  absolute  value  sign,  the  term  becomes  unwieldy 
in  later  expressions  and  is  neglected  at  this  point. 

If  Ip  is  the  moment  of  inertia  of  the  recoiling  mass  about  an 
axis  through  the  mass  center,  we  may  write 

Ir<  $i  +  *  T+  T»  C-13 

In  order  to  evaluate  the  retarding  force,  <^)(t)  .  we  must  know 
what  type  of  resistance  is  furnished  by  the  recoil  system.  The  system  used 
in  the  XM70  launchers  contains  a  preloaded  recuperator  or  air  spring  to 
hold  the  mass  in  battery  position  at  angles  of  elevation.  Sufficient  accuracy 
will  be  obtained  if  a  linear  spring  rate  is  used  to  represent  the  recuperator. 
The  recoil  force  is  furnished  by  forcing  oil  through  a  variable -area  orifice. 
The  total  retarding  force  is  then  made  up  of  the  following  terms:* 

K,  k  (u  -  u),  R  (u^  ,  u),  H.,  H,  (u^  ,  u)  _  C-14 

o  i  i 

The  nomenclature  used  is 


recuperator  preload 

recuperator  spring  rate 

hydraulic  resistance 

constant  friction  force 

hydraulic  system  seal  frictional  force 

in-fcattery  position  of  the  recoiling  mass 


*  For  a  more  complete  description  of  the  recoil  system,  see  Bimonthly 
Report  No.  7  of  this  Project  (8130) 
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Fig.  C  -5  FREE  BODY  DIAGRAM  OF  THE  CRADLE 


The  cam  reaction  force  In  the  direction  of  recoil  may  be  represented  by 

C  (ii,  u,  u);  the  actual  functions  will  be  introduced  later,  in  order  to  keep 

the  lengths  of  the  expressions  to  a  minimum.  The  total  retarding  force  is,  then 

*2  .  .  ..  ...  , ».2 


<£ 


(t)  =  K  +  k  (uq  -  u)  +  eR  (u  ,  u)  +  eH^  +  eH^  (u  ,  u)  +  efC  (u,  u,  u). 

C-15 

In  this  expression,  f  =  0,  except  when  u^u-^u^  and  u<0,  that  is  when 
the  recoiling  mass  is  counterrecoiling  through  the  cam,  beginning  at  Uj  and 
ending  at  u^,  otherwise  f  a  +1, 

We  will  now  consider  a  rigid  cradle  with  a  free  body  diagram  as 
shown  in  Fig.  C-5.  It  is  assumed  that  there  is  no  relative  angular  motion 
between  the  cradle  and  recoiling  mass;  because  of  this,  the  coordinates  of 
the  mass  center  of  the  cradle  may  be  immediately  written  from  Eq.  C-l. 

This  is  done  by  substituting  the  constant  length  for  u.  These  equations 


are 


and 


xc  =  JLcob  (  ft*  0,)  +  bj  cos  G 
yc  =  X  sin  {  J3+  0^)  +  bj  sin  L 


C-16 


The  length  bj  is  thb  distance  along  the  u-axis  from  the  trunnion  to  the 

cradle  mass  center,  and  x  and  y  are  the  absolute  coordinates  of  the 

c  'c 

mass  center  of  the  cradle.  The  corresponding  velocity  and  acceleration 
components  may  be  found  by  differentiation  of  Eq.  C-16  or  by  similar 
substitution  of  b^  for  u  in  Eq.  C-2,  3,  and  4.  The  components  of  acceler¬ 
ation  of  the  cradle  mass  center  are: 


xc  =  -  /0l  sin(/J+01)  -  JLxj)*  cos  (/}+  0t) 

-  b,  (  <t\  -  $2>  sin  Uj  -b  1(01  -f  02^2  cos  (jJ j 

Y c  ~  JL  cos  ( J3+  0L)  -l<pz  sin  (J3+  0l> 

cos  -b,<0i  +  02>  sin 


C-17 


C-1S 


If  we  let  the  vector  sum  of  the  comoonents  of  x  and  y  be  resolved  into 
v  *  c  c 

components  ^  and  .  respectively,  in  the  u-,  v-coordinate  system  using 

Fig.  C-3,  we  have 

£  =  -  £  sin  -  X0  cos  U7 2  -  bj  (  0  J  +  0?)  C-19 
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and 

X  s  '-L  COB  {Jz~  L  sin  Uz  -  bj  (  (f)  J  +  Cp2)  C-20 

The  equation*  of.  motion  of  the  cradle  mas*  center  may  now  be  fotmd  by 
summing  force*  in  the  u-  and  v-directions.  These  are: 

mc  £*  =  lj/(t)  -  (|)(t)  -  S  coa  J2  -  mc£  sin  U j  C-21 

mc  X  ~  V  (t)  -  N  -  mcg  cos  Ui  m  s  sin  f 2  ,  C-22 

where  the  following  nomenclature  is  used: 

(t)  trunnion  reaction  in  the  u-direction 
V  (t)  trunnion  reaction  in  the  v-direction 
S  force  in  the  elevating  system. 

The  displacement  of  one  end  of  the  elevating  system  relative  to 
the  other  is  given  by  b(j£)2  sin  f2.  The  force  developed,  assuming  a  linear 
spring  and  damping  proportional  to  the  first  derivative  of  the  displacement,  is 

S  =  K2  +  k2b(sin  T2)02  +  c2b(sin  ?^)02  C-23 

The  preload  is  given  by  K2#  the  rate  by  k,  and  the  damping  coefficient  by 
c2-  In  this  expression,  the  angle  f 2  1®  assumed  constant.  The  angular 
equation  of  motion  of  the  cradle  is  found  by  summing  moments  about  the  mass 
center.  From  this  we  get 

Ic(<^l+<?)2)  =  -bIV^-T  +  N  (bi  -  «)  -  S  (b  -  bj)  sin  J2  C-24 

In  Eq.  C-24,  I£  is  the  moment  of  inertia  of  the  cradle  about  the  mass  center. 
The  preload  K2  is  the  preload  necessary  to  balance  the  component  of  force 
due  to  the  weight  of  the  recoiling  mass  and  the  cradle  itself. 

The  final  equation  of  motion  may  now  be  found  by  summing 
moments  on  the  firing  support  about  the  fixed  point  0  (see  Fig.  C-6): 

1Q  CL/j  =  Uf(t)  f,sin  C -  V  (t)  J^cos  U ^  ~  S  A,  sin  'f ^ 

‘  K1  "  kl  01  "  C1  01  C-25 

For  any  angle  of  elevation,  it  will  be  assumed  that  f  3  is  a  constant;  also 
the  following  relationship  will  be  used: 
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7  3  ~  'h  ‘0(2  J0  -  (X 


C-26 


a=o 

We  must  now  express  the  recoil  and  cam  forces  in  terms  of  the  variables  of 
this  analysis,  The  recoil  force  is  given  as 


eR  (i2  ,  u)  - 


111 


2  «  C02<V 


C-27 


where  the  following  nomenclature  is  used: 
density  of  th$  oil 


A 

P 

Co 


piston  area 

orifice  coefficient 

orifice  area  (function  of  u) 


For  this  analysis,  it  is  assumed  that  the  seal  friction  is  directly  proportional 
to  the  recoil  pressure  (or  force)  and  may  be  expressed  as 

C-28 


H, 


6r  . 


where  (5  is  a  constant. 


A  cycloidal  indexing  cam  shape  is  used  in  this  analysis  (as  on  the 


XM70  launcher).  The  axial  component  of  force  is 


I 


B 


c 


u  -  u 


•%  -rr-  U  *  U* 

UJL.  si„  zTC - 1 


u . 


u 


U2  ’“I 


(  1  -  cos  ZTT- — 7-—-)  +  (!  -  cos  2  7T 


u  -  u 


w 


C-29 


‘2  “1  “2  “1 
The  launcher  has  two  breech  clusters  containing  three  rounds  each.  For  this 
analysis,  it  is  assumed  that  only  one  cluster  exists  containing  six  rounds; 
however,  the  proper  moment  of  inertia,  Ijy  corresponding  to  the  combined, 
double  cluster  will  be  used.  In  Eq.  C-29.  h  is  the  height  of  the  cam  path 
development  and  is  the  pitch  radius. 

Substituting  the  last  three  relationships  into  Eq.  C-il,  12,  13, 

21,  22,  24,  and  25,  and  eliminating  the  four  unknown  reactions,  T,  V,  N, 
and  lj/,  leaves  three  equations.  These  three  equations  are  second  order, 
non-linear,  differential  equations,  with  the  variables  (^)j.  (^D^,  anc*  u' 
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■/.  J- 


I! 

I! 


f! 


|  i! 

I ? 

r 

M 

■  y  i  * 

:  li 


I 


They  are: 

mr  (  /  »in  U2  +  u)  <^j  + 


ef  — - 


R1  (u2  *  ul> 


2  (1  -  cos  1% 


u  -  u 


u2  “  U1 


1  .2 

— ;  -m 


-mru  (  <$x  +  (^)2)2  +  F  (t)  -  K>  k  (  ufl  -  u)  +  mtl  cos  U2  fy 


-  e 


u  -  u. 


(»>($)(-?!_ )  t/ia,  (M>2-,  ,ln , 


(1  -  cos  2  7T 


2«  Co  «o 

— I-  ) 
U2  "  U1 


and 


<  “2  -  “!> 


u2  "  U1 


•2 


u  -eHj  +  mrg  sin  U j 


C-30 


01  £  *o  +  mc^  +  l£coa  ^2  +  rnr^COB  +  mc  kj/cos  U2J  + 
02  |^mr  u/cos  U2  +  mc  b^cos  CJ2j  +  u  ^-ef/pin  (J2* 

7T - )2  1  =  J  ain  LA  -  b.  m  • 

U2-U1  J  2L  1  c 


h 

- 5 - — — t"  { 1  -  cos  2 

R1  <u2  “  “l> 


•(  <^)j  +  0,)  +  m^cos  U2  +  K  +  k  (uq  -  u)  +  mcg  sin  L/  j  + 
/cos  U2  [-  2  mru  ( <$)  j  +  02)  -  g  cos  (Jj  (mr  +  mc)  | 

r  -  A  2  C 

+  e  /  sin  LA  ( 1  +  (5  )  - i---E — - 

«-  2g  c/0(  2 


+ 

2 


u  -  u 


o  o 

U  -  l<  . 


•  2  JB  2  7T  h‘ 

U  +  H.  +  f  -2_ - _ 

R1  (u2  -  V 3 


sin  2  7T  —  1  (1  -  cos  277" - —  )  u 


U-»  -  u 
2  1 


u2  ‘  U1 


-ci  0i 


and 


~kl  01  -  K1  +  [K2  +  k2»(sin  T2)02  +  C2  b(sin  72)0. 

•  /in  (  U2  -  J  2)  -  /'  sin  (  70  -  C< ) 

0!  Ic  +  Ir  r  mr  u  {/cos  LA  +  u)  +  bt  mc  (/cos  (J2  +  bj 
+  02 


C-31 


I  +I  +  m  u  +m  b,“ 
c  i  r  cl 


-  0j2  /3i«  L/2  (mr  u  +  mc  bj) 
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.  [ 


i. 


hli 

!i! 


I 

I 


.  r.  <y 

H 

I  T 


y 

<*% 

I 

i 

l 


2  mr  u  u  (  0^  +  (^)2)  -  g  cos  Ul  (mr  u  +  mc  bj)  +  C*F(t) 


b  sin 


72  [  k2  +  b  .In  /2  <K2  02  *  c2  4>2)1 


C-32 


III.  DISCUSSION 

Bated  upon  the  assumption  maae,  tne  solution  to  the  three 
differential  equations  of  motion  just  derived  will  furnish  the  first  mode 
response  of  the  launcher  in  the  vertical  plane.  The  reaction  forces  and 
moments  may  also  be  found  when  the  equations  of  motion  are  solved.  Because 
of  the  complexity  of  these  expressions,  their  solution  is  obtained  on  the 
Technology  Center  Univac  1105  digital  computer.  The  programming  of  these 
equations  and  the  results  obtained  is  presented  in  other  appendices. 


tlnOUD  USCUCU  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


C-  1 » 


APPENDIX  D 

DIGITAL  COMPUTER  PROGRAM, 

THREE-DEGREE -OF-FREEDOM  MATHEMATICAL  MODEL 

R,  M.  Brach 

I.  INTRODUCTION 

The  equations  of  motion  of  a^  three  -degree  -of-freedom  system 
representing  an  XM70-type  launcher  are  presented  in  Appendix  C.  Because 
of  the  complexity  of  these  equations,  they  are  being  solved  numerically  on 
Armour  Research  Foundation's  Univac  1105  digital  computer.  In  this 
appendix,  the  computer  program  Is  described  and  some  results  are  compared 
with  previous  work. 

II.  COMPUTER  PROGRAM 

Basically,  the  computer  program  constitutes  the  numerical 
solution  of  a  set  of  three,  second-order,  non-linear  differential  equations. 

The  solution  was  programmed  in  a  non-algebraic,  symbolic  language 
called  USE  (Univac  Scientific  Exchange).  A  side-by-side  listing  of  the 
program,  as  written,  and  its  octal  representation  are  shown  in  Fig.  D-3. 
Program  output  format  is  illustrated  in  Fig.  D-4,  which  is  the  first  page 
of  output  of  a  typical  run. 

A  subroutine  was  wri  tten  for  general  use  by  Robert  Floyd  of 
Armour  Research  Foundation  based  on  the  numerical  technique  of  Runge  and 
Kutta  for  the  solution  of  up  to  50,  first-order,  simultaneous  differential 
equations.  It  is  applied  here  to  solve  the  three,  second-order  equations  after 
a  reduction  of  their  order  by  three  triviai  substitutions,  such  as  u  =  s.  The 
subroutine  is  shown  listed  after  the  main  program  with  the  call  address, 
RWF999.  One  other  special  subroutine  was  written,  RWF998,  to  allow  the 
use  of  two  tabular  functions:  in  this  case,  an  experimental  powder -gas 
forcing  function  and  a  table  of  recoil-system  orifice  areas.  This  subroutine 
employs  linear  interpolation  between  points.  It  is  also  listed  after  the  main 
program.  All  other  subroutines  used  are  from  the  standard  ARF  library. 
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D-l 


Fig.  D-l  shows  a  gene  ral  program  flow  diagram.  A  detailed 
flow  diagram  showing  how  the  time  of  application  of  the  powder  gas  force 
is  derived  from  the  position  and  velocity  of  the  recoiling  parts  is  shown  in 
Fig.  D-2. 

III.  RESULTS 

Two  comparisons  of  the  output  of  this  program  will  be  made 
in  this  appendix:  0)  a  comparison  with  a  previous  computer  program  of  the 
parent  project  and  (2)  a  comparison  with  experimental  curves  obtained  from 
firings  of  the  Prototype  No.  1  launcher. 

The  physicalUauncher  parameters  used  for  the  first  comparison 
are  taken  from  Prototype  No.  1.  The  response  curve  from  an  early  computer 
program  is  shown  in  Fig.  D-5;  the  corresponding  response  curve  from  the 
present  computer  program  is  shown  in  Fig.  D-6.  Table  0-1  points  out  the 
significant  similarities  of  and  differences  between  the  two  sets  of  equations 
of  motion  and  launcher  parameters. 

The  two  curves  are  seen  to  differ  slightly  in  amplitude  and 
frequency.  The  amplitude  differen  ces  may  be  accounted  for  by  the 
difference  in  recoil  forces  between  the  two:  the  recoil  force  in  the  current 
program  peaks  significantly  near  the  end  of  recoil  as  does  the  recoil  force 
of  the  actual  launcher,  whereas  in  the  early  program,  the  recoil  force 
remains  nearly  constant.  The  frequency  differences  probably  are  due  to  the 
inclusion  of  the  elevating  system  flexibility  in  the  current  program,  even 
though  that  system  is  very  stiff. 

The  second  comparison  utilizes  the  experimental  record  of  the 
first  shot  of  Burst  No.  L68-L73,  which  was  fired  at  Redstone  Arsenal  on 
February  22,  1959  (presented  in  Appendix  II-A  of  Report  No.  1 1  of  the  parent 
project).  Figure  D-7  shows  the  trail  strain,  cradle  motion,  and  recoil  force 
measured  during  that  burst.  The  response  as  given  by  the  current  computer 
program  is  shown  in  Fig.  D-8.  The  input  data  to  the  computer  program  to  be 
compared  are  shown  in  Table  D-2.  The  chosen  elevating  system  stiffness 
value  is  seen  to  be  relatively  low;  it  may  be  recalled  that  the  elevating  screw 
brackets  on  the  Prototype  No.  1  were  found  to  be  relatively  flexible. 
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Table  D-l 

COMPARISON  OF  MODELS  AND  LAUNCHER  PARAMETERS 


Item 

Old  Program 

New  Program 

Degrees  of  freedom 

Recoil  motion,  car- 
raige  rotation,  hop  of 
firing  base 

Recoil  motion,  car- 
ruige  rotation,  cradle 
rotation  relative  to 
the  carriage 

Launcher  Geometry 

Pilot  No.  1  at  15°  . 
elevation 

Pilot  No.  1  at  15° 
elevation 

Value*-  of  Mass 

Pilot  No.  1 

Pilot  No.  1 

Elevating  System 
Stiffness 

Rigid 

l  x  106  lb/ in. 

Equivalent  Carriage 
Stiffness 

100  x  10^  in.  -lb/ radian 

100  x  10^  in.  -lb/ radian 

Recoil  System 

Computed  from  equa¬ 
tions  of  a  theoretical 
recoil  system 

Computed  from  equa¬ 
tions  of  a  theoretical 
recoil  system  using 
actual  orifice  areas 
from  Pilot  No.  3  rods. 

Auxiliary  damping 

None 

None 

Table  D-2 

INPUT  DATA  TO  COMPUTER  PROGRAM 


Item 

Value  in  Program 

Launcher  Geometry 

Pilot  No.  3  at  45°  elevation 

Values  of  Mass 

Pilot  No.  3  (See  parameter  list  in 
Program,  Fig.  D-3) 

Elevating  System  Stiffness 

1.405  x  104  lb/in. 

Elevating  System  Damping  Coeffi¬ 
cient 

90  lb  -sec/in. 

Equivalent  Carriage  Stiffness 

45  x  10°  in.  -Ib/radian 
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The  similarity  between  the  experimental  and  computed  curves 
is  striking,  even  down  to  details  such  as  the  influence  of  indexing  forces 
between  times  0.26  and  0.4.  Especially  important  is  the  reproduction  of 
the  computed  curve  of  cradle  motion  in  the  experimental  curve  Just  prior  to 
the  time  of  firing  of  the  second  round.  This  cradle  motion,  the  most  undesir¬ 
able  part  of  the  gross  response  of  the  Prototype  No.  1  launcher,  was  subse¬ 
quently  reduced  by  designing  a  stiffer  elevating  system.  By  means  of  an 
accurate  computer  program,  however,  the  optimum  combination  of  para- 

o  ' 

meters  can  easily  be  found  for  a  desired  response. 

IV.  CONCLUSIONS 

The  examples  presented  in  this  appendix  indicate:  (1)  that 
the  current  computer  program  corresponds  closely  to  previous  theoretical 
work,  which  means  that  errors  do  not  exist  either  in  the  equations  that  were 
programmed  or  in  the  program  itself;  and  (2)  that  the  data  derived  from  the 
current  computer  program  represent  gross  experimental  launcher  response 
very  accurately,  which  means  that  the  program  may  be  used  to  improve 
that  response. 
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Read 
input  data 


Calculate 

constants 


RWF  999,  Numorica  Integration  Subroutine 

!  I - 1 


Calculate  derivatives 
from  equations  of 
motion 


Calculate  variables 
using  Rungc-Kutta 
procedure 


Print 

output 


Calculate 
powder -gas  force 


Is  time  less  ' 
than  final  time? 


Fig.  D-l  COMPUTER  PROGRAM  FLOW  CHART 
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FIG.  D-3  COMPUTER  PROGRAM  (CONTINUED) 
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FIG.  D-3  COMPUTER  PROGRAM  (CONTINUED) 
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appendix  e 

OPTIMUM  ELEVATING  SYSTEM  PARAMETERS 
Coleen  Murray 


I.  INTRODUCTION 

The  purpose  of  this  study  was  to  determine  the  effects  of 
various  combinations  of  elevating  system  stiffness  and  damping  upon  gross 
launcher  motion  and  to  find  the  optimum  combination  of  such  stiffness  and 
damping.  The  optimum  combination  was  defined  as  that  combination  which 
yields  a  consistently  low  value  of  cradle  displacement  prior  to  the  firing  of 
each  of  the  successive  rounds  of  a  burst  and  which  concurrently  yields  a 
minimum  corresponding  cradle  velocity.  The  work  was  based  on  the  assump¬ 
tion  that  cradle  oscillations  sustained  during  a  burst  firing  contribute  directly 
to  launching  errors  that  cause  dispersion  of  the  shots  of  a  burst.  This  study 
is  concerned  primarily  with  cradle  motion,  although  barrel  motion  relative 
to  the  cradle  is  also  a  major  cause  of  launching  errors. 

This  study  utilizes  the  work  of  Philip  Meyfarth  -1  in  order  to 
define  the  general  area  of  optimum  damping  for  a  given  stiffness.  Optimi¬ 
zation  of  Meyfarth's  third-order  system  yielded  an  optimum  damping  coeffi¬ 
cient  for  any  given  elevating  system  stiffness.  The  damping  determined  by 
Meyfarth’s  work  was  used  to  indicate  a  value  about  which  the  damping  should 
be  varied  in  the  Model  V  solutions  in  order  to  obtain  the  actual  optimum. 

The  ranges  of  optimum  damping  thus  obtained  for  the  several  values  of 
stiffness  were  then  compared  to  each  other,  and  a  range  of  optimum  stiffness, 
with  corresponding  optimum  damping,  was  chosen.  The  third-order  system 
was  used  because  of  the  availability  of  the  response  data  in  the  reference  and 
the  resulting  simplification  of  the  optimization  procedure  for  the  three-degree  - 
of-freedom  computer  solution. 

A  few  variations  of  equivalent  carriage  stiffness  for  the  optimum 
elevating  system  were  also  investigated  to  determine  the  effect  of  this  stiffness 


1/ 


Meytarth.  P.  ,  Dynamic  Response  Plots  and  Design 
Order  Linear  Systems.  Research  Memn  Mr,  R  \a  " 
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on  cradle  motion  and  trail  hop,  which  is  the  upward  deflection  of  the  trails 
observed  during  firing  of  the  XM70E1  launcher  when  the  recoiling  assembly 
is  in  its  extreme  recoil  position.  The  values  of  launcher  parameters  main¬ 
tained  constant  in  the  study  are  those  of  the  XM70E1  Prototype  No.  3 
launcher  at  0°  traverse  and  45°  elevation. 

II.  DISCUSSION 

A.  Third-Order  Analysis 

The  angular  motion  of  the  launcher  was  approximated  by  the 
motion  of  the  following  linear  system  with  the  proper  conversion  of  units: 


VY\ 


Fig.  E-l  MODEL  OF  THIRD-ORDER  SYSTEM 

where 

m  is  an  equivalent  mass  moment  of  inertia  of  the  tipping  parts 

about  the  trunnion,  7500  in.  -lb-sec** 

k  ,  is  the  equivalent  angular  stiffness  of  the  elevating  system 

c  ^  is  the  equivalent  angular  damping  coefficient  of  the  elevating 

system 

k  is  the  combined  angular  stiffness  of  the  carriages  and  trails, 

3.  5  x  10 ~  in.  —lb/ radian 

F(t)  is  a  forcing  function  applied  to  the  trunnion. 
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The  equation  of  motion  for  m  is: 
m  x,  3  -k  2  (Xj  -  X2)  -  c  2  -  *2) 


E- 1 


From  the  statics  at  x^  ,  we  have: 

1  1 

kx2  =  k  2  (Xj  -  x2)  +  c  2  (Xj  -  x2)  +  F(t) 

By  eliminating  X2  from  E.-l  and  rewriting  in  terms  of 
differential  operators,  one  obtains: 

x,  (c  ,)  +  Xj  (k  +  k  j)  +  D  x^  (c  ^  k/m)  +  Xj  (k*k  ^  fm) 


E-2 


1  '2'  l  *  2'  l 

(c'2/m)  D  [f( t)]  +  (k'2/m)  [F{t)] 


E-3 


The  left-hand  side  of  Eq.  E-3  may  be  factored  into  three 

* 

factors,  each  representing  a  root  of  the  system  .  When  two  of  these  roots 
are  complex  conjugates,  the  parameters  ,  (_j,  and  ^  may  be  defined  by 
the  following  equations: 


c  2m 
■ — 1 — 
k  k  - 


ul 


(ktk'2)m  .  (jtru  * ») 


k  k 


1  -2  "■ 


2  U* 

-  lT U  *  2  £  >  • 

U 


E-4a 


E-4b 


E-4c 


where 


T  is  the  constant  of  the  factor  representing  the  real  root 

£  is  the  damping  ratio  associated  with  the  complex-conjugate  pair 

U  is  the  undamped  natural  frequency  associated  with  the  complex- 
conjugate  pair. 

Let  the  right-hand  side  of  Eq.  E-3  be  an  impulse  function  of  strength  I; 
then,  the  solution  of  Eq,  E-3  is: 

For  a  complete  discussion  of  the  following  technique  consult  Meyfarth, 
Reference  1 . 
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The  first  of  the  two  terms  in  the  right-hand  side  of  Eq.  E-5 
is  an  inverse  exponential  which  decreases  with  increasing  time;  the  second 
is  the  product  of  an  inverse  exponential  and  a  sinusoidal  function.  The 
sinusoidal  function  will  decay  in  the  shortest  time  when  the  values  of  the 
coefficients  of  t  in  the  exponents  are  minima.  The  values  of  T,  LJ ,  and 
£  which  appeared  to  yield  the  best  damping  were  chosen  from  the  graphs  in 

Reference  1  and  the  corresponding  value  of  damping  was  then  evaluated. 

.1 

Figure  E-2  shows  the  response  for  several  values  of  k  ^  with 
the  corresponding  optimum  damping.  The  angular  stiffness  and  angular 

i  i 

damping  coefficient,  k  ^  and  c^,  were  converted  to  values  corresponding 
to  linear  restraints  by  multiplying  by  10“^  (the  inverse  of  the  square  of 
the  distance  from  the  trunnion  to  the  elevating  screw  bracket). 

The  following  comments  are  made  on  the  statement  proceeding 
Equation  E-4a,  concerning  the  pair  of  complex  roots.  Large  values  of 
stiffness  always  result  in  two  of  the  roots  being  complex  conjugates.  However 
below  a  certain  stiffness  it  is  possible  to  get  three  real  roots  for  some  ranges 
of  damping.  Furthermore,  as  the  stiffness  decreases,  there  is  an  increase 
in  the  range  of  damping  which  causes  this.  When  these  damping  values  were 
used  as  input  to  the  computer  solutions,  they  resulted  in  a  continuously 
increasing  (non-oscillatory)  displacement,  as  shown  in  Figure  E-3. 

B.  Computer  Solutions 

Initially,  the  computer  solutions  were  investigated  for  the 
combined  carriage  and  trail  stiffness  of  45  x  10 ^  in.  -lb/radian,  which 

is  about  30%  higher  than  the  measured  stiffness.  A  range  of  vary  good 

4  ... 

response  was  found  for  stiffness  of  1.0  x  10  Ib/in.  and  damping  of  3o0  to 
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540  lb-sec/in.  Several  graphs  of  cradle  rotations  resulting  from  paired 
values  of  stiffness  and  damping  in  this  range  are  shown  in  Fig.  E-4.  The 
choice  of  an  optimum  combination  here  would  be  arbitrary  because,  at  the 
time  of  each  succeeding  firing,  the  cradle  displacement  is  uniform  (Fig.  E-5) 
and  the  velocity  is  small  for  all  of  the  responses, 

To  illustrate  that  this  optimum  system  is  an  improvement,  its 
response  is  compared  in  Fig,  E-6  with  a  response  similar  to  that  of  the 
present  weapon  (with  little  damping).  Also  shown  is  an  over -damped  response, 
to  demonstrate  that  an  optimum  damping  does  exist. 

It  was  realized  that  there  were  limits  on  the  obtainable  values 
of  stiffness  and  damping  in  the  elevating  system.  Yery  high  stiffness  was 
shown  to  be  difficult  to  achieve  in  previously  built  launchers,  and  extremely 
low  stiffness  presents  a  problem  of  stability  during  sighting  and  also  is 
subject  to  excessive  preload  deflections.  Comparison  of  experimental  and 
analytical  results  indicates  that  damping  corresponding  to  at  least  90  lb-sec/in. 
is  inherent  in  the  launcher  system;  thus,  no  smaller  damping  is  feasible. 

Next,  a  much  lower  stiffness  with  increased  damping  was 
investigated  to  discover  if  cradle  rotation  would  thus  be  decreased.  Fig.  E-3 
shows  that  for  damping  between  the  non-oscillating  range  up  to  the  maximum 
damping  used,  increased  damping  gives  better  response.  Nevertheless, 
these  responses  are  not  as  good  as  those  with  the  higher  spring  rate  (see 
Fig.  E-7  for  comparison).  Also,  the  study  of  a  4-rd  burst  (Fig.  E-8)  showed 
that  cradle  rotation  increased  with  each  round. 

Finally,  the  computer  solutions  were  investigated  for  the  com¬ 
bined  carriage  and  trail  stiffness  of  35  x  10  ^  in.  -lb/ radian;  this  is  very 
close  to  the  experimental  value  for  Prototype  No.  3.  Although  the  initial 
maximum  cradle  rotation  was  somewhat  greater  than  for  the  stiffer  structure, 

the  range  of  good  response  was  similar.  The  optimum  response  still  occurred 

4  a 

for  a  stiffness  of  1.0  x  10  to  1.  4  x  10  lb/in.  and  damping  of  360  to  540 
lb-sec/in. 

Trail  hop  was  significantly  affected  by  the  carriage  stiffness, 
kj,  as  shown  in  Fig.  £-9,  which  is  a  plot  of  the  cradle  rotation  for  three 
systems  with  otherwise  identical  parameters.  The  maximum  negative 
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rotation  at  about  0.  2  sec  of  the  supporting  structure  is  an  indication  of  the 
relative  magnitude  of  the  hop,  illustrating  that  trail  hop  greatly  decreases 
with  increasing  carriage  stiffness. 

During  the  course  of  the  study,  it  was  noted  that  the  optimum 

damping  determined  from  the  computer  solutions  agreed  with  the  optimum 

.4 

predicted  by  the  third-order  analysis,  where  linear  stiffnesses  from  1.0  x  10 
to  1.4  x  10^  lb/in.  corresponded  to  optimum  damping  from  460  to  510  lb-sec/in. 
Such  close  agreement  was  interesting  and  unexpected,  in  that  the  Model  V 
equations  include  the  actual  input  force,  cam  path,  and  many  nonlinear  launcher 
effects  not  considered  in  the  third-order  system. 

111.  CONCLUSIONS 

On  the  basis  of  the  investigations  described  above,  it  is  con¬ 
cluded  that  improved  firing  accuracy  can  be  achieved  for  the  system  studied 

by  adding  damping  to  the  elevating  system.  Specifically,  the  stiffness  of 

4  4 

the  elevating  system  should  be  about  1.0  x  10  to  1.4  x  10  lb/in.  and 
the  damping  of  the  elevating  system  should  be  about  360  to  540  lb-sec/in. , 
for  the  maximum  improvement. 
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Fig.  E-2  CRADLE  ROTATION  FROM  THIRD-ORDER  SYSTEM  ANALYSIS 


Fig.  E-5  CRADLE  ROTATION  FOR  A  6-RD  BURST  WITH  OPTIMUM  ELEVATING  SYSTEM 


CRADLE  ROTATION  FOR  A  4-RD  BURST 
WITH  LOW  ELEVATING  SYSTEM  STIFFNESS 


APPENDIX  F 

GENERAL  EQUATIONS  OF  MOTION  OF  LAUNCHER  DYNAMICS 


I.  INTRODUCTION 

The  3-degree -of-freedom  model,  presented  in  Appendix  C, 
was  programmed  for  a  digital  computer  and  has  been  used  extensively. 

As  a  3-degree -of-freedom  model,  it  has  the  advantage  of  relative  simpli¬ 
city  and  also  reproduces  much  of  the  primary  behavior  of  tho  XM70-typo 
launcher.  On  the  other  hand,  the  very  simplicity  precludes  its  application 
to  many  of  the  presumedly  secondary  features  of  the  launcher.  In  order  to 
study  these  features,  a  more  complete  model  with  up  to  n-degrees  of 
freedom  was  constructed. 

The  following  are  some  of  the  salient  features  of  this  more 
complete  model  that  are  present  in  addition  to  the  features  of  the  3-degrce- 
of-freedom  model: 

1.  The  base  of  the  launcher  is  not  assumed  to  be  fixed  to  a  rigid 
foundation;  instead,  the  trail  pads  and  base  plate  are  placed 
on  visco-elastic  supports  which  can  be  assigned  parameter 
values  corresponding  to  various  terrain  compositions  (clay, 
sand,  concrete,  etc.).  The  model  takes  cognizance  of  the 
fact  that  this  may  be  a  one-way  support  condition  so  that  trail 
hop  and  base  hop  are  directly  reproduced.  Should  the  need 
arise,  the  model  can  be  easily  altered  to  reproduce  progressive 
trail  dig -in  and  base  settlement. 

2.  It  is  recognized  that  the  launcher  does  not  recoil  straight  back, 
maintaining  its  position  in  its  initial  vertical  plane.  The 
tendency  for  three-dimensional  rotation  about  the  ball  joint  is 
allowed  for.  Reaction  forces  and  moments  at  the  ball  joint, 
the  trunnion  supports,  the  cradle,  and  the  recoiling  parts  are 
admitted  as  three-dimensional  quantities. 

3.  The  flexibility  of  the  recoiling  assembly  is  provided  for,  as  are 
the  lateral  and  vertical  variations  of  center  of  gravity  and  iner¬ 
tial  force  location  from  shot  to  shot.  The  additional  features 
allow  the  reproduction  and  study  of  the  ''left-right  sequencing", 
phenomenon  observed  in  tests  of  the  XM70  prototype  launchers. 

4.  The  transmission  of  rifling  torque  to  the  launcher  structure  is 
included,  as  are  the  force  parameters  and  geometrical  freedoms 
necessary  for  interpretation  of  its  effect. 
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5.  The  rifled  firing  tube  is  represented  as  a  flexible  piece  with 
many  degrees  of  freedom  in  both  the  vertical  and  horizontal 
planes.  This  permits  the  investigation  of  tube  whip  and 
vibrations  and  their  effect  on  both  the  supporting  structure  and 
the  projectile. 

6.  The  projectile  itself  is  a  part  of  the  model.  The  resultant 
effects  include  the  variation  of  mass  during  firing  and  the 
influence  of  traveling  forces  on  launcher  tube  motion. 

7.  A  result  of  the  inclusion  of  the  increased  geometrical  freedom 
and  applied  forces  is  the  introduction  of  Coriolis  forces  (in 
three  dimensions)  and  gyroscopic  moments  not  previously 
considered.  An  additional  result  is  the  inclusion  of  more  of 
the  coupled  effects  of  the  motion  of  the  various  launcher 
components. 

This  appendix  presents  a  summary  of  the  derivation  of  the 
basic  equations  of  this  model,  as  well  as  graphical  and  verbal  definitions 
of  the  known  and  unknown  parameters  involved  with  the  basic  equations 
that  govern  the  model.  The  motion  of  the  launcher  can  be  best  expressed 
in  terms  of  a  set  of  interdependent  right-hand  three-dimensional  coordinate 
systems,  as  shown  in  Fig.  F-l.  Coordinate  system  1  (  [s^]  )  is  fixed 
in  the  earth  at  the  trail  pads  with  y^  vertical.  Then,  Xj  is  defined  by  a 
plane  through  y^  and  the  initial  tube  centerline.  Systcrps  t  S^l  and 
have  origins  at  the  ball  joint  and  arise  from  a  translation  and  a 
rotation  through  angle  (0  j  +  0  .,)  and  the  x^,  y^  plane.  Angle  &  ^  is  a 
measure  of  trail  flexibility,  while  angle  0  ^  os  a  ball-joint  rotation. 
System  \_S^]  arises  through  a  rotation  through  angle  0^  about  y^.  Thus, 
y^  coincides  with  y^.  System  ^Sg]  arises  through  a  rotation  through 
angle  0^  about  x^;  Xj.  coincides  with  x^.  It  is  apparent  that  &  y  0 
0  j,  represent  freedoms  at  the  ball  joint.  Wearrive  at  [  S^j  by  means 
of  a  translation  through  distance  H  ^  and  a  rotation  through  angle 


axis  parallel  to  .  In  addition  to  coordinate  systems.  Fig.  F-l 

illustrates  many  of  the  quantities  which  are  defined  or  following  pages 
and  which  appear  in  the  equations. 
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C*.  +  0  ^  ;  the  origin  of  LS^j  is  at  the  trunnions.  The  angle  is  the  angle 
of  elevation  and  0  is  the  freedom  in  elevation  permitted  by  the  flex 
system.  The  tube  base  is  located  a  distance  J.  along  the  x,  axis.  Left- 
right  sequencing  forces  the  tube  to  swing  through  an  angle  0  ^  about  an 


*mm*<*y 


II.  SPECIFIED  QUANTITIES 


trail  length 
carriage  side  length 


trunnion-to-cradle  center  of  gravity 

ball  joint  to  flex  attachment 

y^-position  of  breech  center  of  gravity 

z^-position  of  breech  center  of  gravity 

distance  from  brooch  center  of  gravity  to  tube  base 

tube  length 

initial  x^-position  of  breech  center  of  gravity 
elevation  angle 

initial  projectile  position  with  reference  to  tube  base 
distance  from  trunnion  to  flex  attachment  point 
horizontal  (z^)  distance  of  cluster  axle  from  tube  centerline 
carriage  side  orientation 


flex  system  orientation 
tube  moment  of  inertia  =  ^ 


TT  .  4  4. 

T  rz  '  ri  * 


tube  mass  density  per  unit  length 
projectile  mass 


J  =  J 

y  2 


mass  polar  inertial  moments  of  projectile 


(0-  A 


prescribed  projectile  spin  history 


breech  mass 


ARMOUR  RESEARCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


JBx 

JBy 

JBz 
m  * 


I 


Ax 


I 


Ay 


Az 

m 


1 


]xc 
\c 
Xzc 
mB 

m 

r 

kr 

k2’  ^  2 


5 

k6’  *  6 
E 

h 

F{  t) 

k8 

V 


mass  polar  inertial  moments  of  breech 

cradle  mass 

mass  polar  inertial  momonts  of  cradle 

concentrated  mass  at  rear  of  trails  (one -half  trail  mass) 

mass  polar  inertial  moments  of  carriage  about  ball  joint 

one-half  mass  of  trails  plus  mass  of  carriage  and  trunnion 
sides  concentrated  at  ball  joint 

mass  of  right  (-f  zq)  cluster  and  projectiles 

mass  of  left{~  z^)  cluster  and  projectiles 

spring  and  damping  constants  at  trail  rear 

spring  and  damping  constants  at  trail  front 

angular  trail  and  carriage  stiffness  (z-axis) 

carriage  stiffness  for  torque  about  y-axis 

carriage  stiffness  for  torque  about  x-axis 


spring  and  damping  constants  for  flex  system 

7 

Young's  modulus  (3  x  10  ) 

coefficient  of  friction  between  breech  and  cradle 
resultant  "powder-gas  force" 
breech  "stiffness"  for  right  -left  sequencing 
eigen  values 
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.."V-  •  •'  V,  .  i’-V  ■ 

V  I 

b 

n 

4  P  2 

wave  numbers  b  =  -Arr—  V 

n  El  n 

I 

1 

:gt\ 

L  bj  =  1.875 

L  b2  =  4.694 

L  b3  ss  7.855 

Lb,  -  10.996 

L  bg  =  14.  137 

|  - 

L  b6  s  17.279 

K 

recuperator  preload 

. .  » 

;  :|.l 

k7 

recuperator  spring  rate 

constant  recuperator  friction  force 

1 

r 

1  i 

oil  density 

piston  area 

C 

o 

| 

orifice  coefficient 

- 

OL 

orifice  area 

• 

g 

acceleration  of  gravity 

>  b 

proportionality  constant;  seal  friction  to  recoil  force 

1  lB 

cluster  inertia  moment 

R 

,  P 

cam  pitch  radius 

!  v 

height  of  cam  path  development 

i  Ul 

initial  cam  engagement  point  ^ 

\  ui  -  u  -  u2 

final  cam  engagement  point  j 

1  “2 

1  III.  LOGICAL  SYMBOLS  * 

j  o  for  4  0 

H{  ^  ) 

-  j  (  \  generic) 

I  1  for  0 

i  for  u  <  0 

£ 

C 

=  '  -1  for  u  >  0 

\  0  for  u  =  0 
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I  At 

IV. 

yT 

yB 

u 

X, 


S 

V 

w 

h 

®3 

^  4 

V 

hn(t) 

7 

«8 

R 


absolute  value  of  A  (A  generic) 

PHYSICAL  VARIABLES  (Unknowns) 
trail  end  deflection 
ball  joint  deflection 

breech  displacement  along  cradle  measured  from  trunnion 
displacement  of  tube  base:  »  u  + 

position  of  projectile  measured  from  tube  base;  positive 


=  V 

y 


6x 


R 


6y 


R/ 


)Z  ^ 


*6x 


'6y 


! 

j 


in  direction 
y^ -displacement  of  tube  at  point  x 
t,  -displacement  of  tube  at  point  x 

angle  due  to  "rigid"  base  movement  Q  ^  *  arc  sin  - ■£- — 

flexural  rotation  of  combined  carriage  flexibility 

transverse  angle  of  rotation  about  y^ 

angle  of  rotation  about  x^ 

elevating  system  flexibility 

right-left  sequencing  flexibility  of  breech 

time  history  for  nth  vibrational  mode  of  tube  in  y^ -direction 

time  history  for  nth.  vibrational  mode  of  tube  in  z,. -direction 

*!  "3  w 

•rotation  about  z„g  .  .  .  ip 


rotation  about 


/3  •  ••  •  W 


3  x 
d  v 
d  .x 


ball  joint  force  reactions  at  trunnion;  positive  acting 
on  carriage;  components  in  system  6 


flex  system  force  components;  positive  acting  on 
carriage;  components  in  system  5 
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trunnion  moment  reactions;  components  in  system  6; 
positive  acting  on  carriage 


force  reaction  between  cradle  and  breech;  components  in 
system  6;  positive  acting  on  cradle 


N 


cx 


cy 


N 

c 

NCZ. 

sx(o,t)  ) 

Sy(0,t) 
My(0,tn 
M  (0,t) 

mm 

N  'l 

py 

N 

pz 
M 

px 
M 

py 

M 

pz 


px 


moment  reaction  between  cradle  and  breech;  components 
in  system  6;  positive  acting  on  cradle 


shear  force  between  breech  and  tube 


moments  between  breech  and  tube 


force  reaction  between  projectile  and  tube;  components 
in  system  6;  positive  when  acting  on  projectile 

moment  reaction  between  projectile  and  breech,  when 
projectile  is  in  breech;  components  in  system  6;  positive 
when  acting  on  projectile 

moment  reaction  between  projectile  and  tube  when  projectile 
is  in  tube  --  transferred  rigidly  to  tube -breech  connection;, 
component  in  system  6;  positive  when  acting  on  projectile 

cos  ( t>C  + 
sin  { 

COS  ip  ' 

sin 

ground  reaction  and  gravity  force  at  y^, 
ground  roack,jn  and  gravity  force  at  y^ 
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V.  DERIVATION  OF  EQUATIONS  OF  MOTION 

The  equations  of  motion  were  derived  with  the  use  of  standard 
matrix  techniques.  These  procedures  are  fully  covered  in  ’’Elementary 
Matrices  and  Some  Applications  to  Dynamics  and  Differential  Equations" 
by  R.  A.  Frazer,  W,  J,  Duncan,  and  A,  R.  Collar,  Cambridge  University 
Press,  1957. 

It  is  well  known  that  the  motion  of  a  body  can  be  represented 
by  equations  which  characterize  the  motion  of  the  center  of  gravity  of 
the  body  together  with  those  which  characterize  the  rotation  of  the  body 
about  its  center  of  gravity.  Figure  F-2  illustrates  tho  general  case.  The 
coordinates  X,  Y,  Z  are  fixed  in  space,  while  x,  y,  z  are  fixed  in  the 
body.  The  motion  of  point  P  is  defined  through  the  rate  of  change  of  R 
or,  equivalently,  the  rates  of  change  of  ^  and  7  .  If  the  point  P  is 
fixed  in  the  body,  G,  the  rate  of  change  of  r  is  completely  defined  by 
the  rotation  of  the  body. 

R  is  normally  referred  to  X,  Y,  Z,  as  is  ^  .  On  the  other 
hand,  r  is  normally  expressed  in  terms  of  x,  y,  z.  If  A  is  the  matrix 
for  transformation  from  X,  Y,  Z  to  x,  y,  z  wo  then  have 

r  =  A  (R  -j*) 

with  respect  to  the  moving  coordinates  the  velocity  of  P  is 
A  R  =  A  J  +  7  +  wr. 

where  w  is  the  rotational  velocity  matrix,  defined  through 
w  =  —  A  A  * 
where  A  *  is  the  inverse  of  A 


R  = 


It  should  be  noted  that 

dR 

at 


is  the  velocity  of  P  with  respect  to  the  fixed  coordinates  and  that  A  R 
merely  is  the  result  of  changing  the  component  bases. 

Differentiation  leads  to  the  following  expression  for  the 
acceleration  referred  to  the  moving  axes 
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i  _  _  1  _ 

A  R  =  A  J  f  2  w  7  +  w  r  +  w  r  +  F 

The  kinematic  relation*  can  be  used  directly  to  obtain  the 
equations  of  motion  for  a  point  mass  or  the  center  of  gravity  of  a  rigid 
body.  In  the  moving  coordinate  system 

MAR  > T 

where  Mis  the  mass  and  F  is  the  force  vector.  For  the  rotational 
equations  of  motion  we  observe  that  the  angular  velocity  can  be  represented 


w 

or 

0 

-w 

2 

w_ 

0 

.  z 

•w 

w 

y 

X 

r“*  3  {w  ,  w  ,  w  ) 
x  y  * 

Let  T  be  the  polar  moment  of  inertia  matrix  of  the  body  represented  by 
the  body  axes.  Then  the  angular  momentum  vector  in  the  moving  system  is 

R  =  7  p 

and  the  time  rate  of  change  of  angular  momentum  is 

mm  *  _  '  — 

HsflsJp+Jp+wJp 
where  M  is  the  vector  of  applied  moments. 

Coordinate  Transformations 

The  applicable  coordinate  systems  are  shown  in  Fig.  F-l. 
These  systems  are  related  through  the  following  transformations. 

^2  =  -1  "  T21 


^  21  =  *  ^  r  yB' 


^3  =  A32  ^2 


cos  (  § 

J. 

1  ' 

^  Z1 

.  n 

sin  (  v  ^  + 

9  ,) 

0 

-sin  (  ^ 

1  + 

*2) 

cos  (  $  J  + 

Q 

0 

0 

0 

1 
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F-.l.O 


R4  ’  *43  *3 


*43 

s  - 

COS  $  a 

0 

-sin 

«3 

0 

1 

0 

sin  $  3 

0 

cos 

s3 

R5  *  *54 

R4 

*54 

n 

1 

0 

0 

fl4 

0 

cos  84 

sin 

0 

-sin  &  A 
* 

cos 

*4 

R6  *  *65  |R5  '  T65( 

765  *  l  t-  2  c0*  Pr  4,ln  Pv  °] 


*65  * 

COS  (  d  +  d  5> 

sin  (  C<  + 

85)  o 

-sin  (  0 (  +  ^ 

cos  (  o(  + 

d  5)  o 

0 

0 

1 

R7  ‘  R6  -  *76 

5  (Xj  *  ^  »  0) 

R8‘  =  R7  "  ^87 

Ta?  =  (x  .  0  .  0) 


A  K 
98  8 


r 


98 


cos 


V, 


-sin  ^  cos  ^  0  . 


0  0  l 


ARMOUR  RESEARCH  POUNOMION  Of  ILLINOIS  .INSTITUTE  Of  TECHNOLOGY 


t  n-yjir  w*  .  jl  -;f ' *f!  •  »£,>  j  ^ 


F-  12 


V1,  AUXILIARY  quantities  and  expressions 

Aw  s  14,  (o)  =  2.  / ( p  L 

a  V‘“(cA=-  •=  ^  UL-t-^akiL  \ 
D”  ^  f  <7U*  \s»*M!cUj 

(f )  =  =g“  lo^-cos,  + 

Vnt:)5  |=-^,Mrt  U, c  +,s>-o  t=nl  *■ 


.  Kosg  UL»o»LLV  , 

\  V*4  L+i  M  L„  “ 


g 


ao 


<wflvz  y,u)^t) 

~ii.i)  =  i  v»(i)u.,d) 


w*  I 

OO 


^(t) *  io-;4  vAtv^ct' 

,r^ = |  t-l  f 'u-)k(t) 

^  - 1  * 


MOUS  «  E  S  E  A  »  C  M  FOUNDATION  OF.  lUlNO 


15  INSTITUTE  OF 


*M] 


TECHNOLOGY 


H  •  I 

6»M  s  Vtin)  =  o/(i)*|  T„'(j)q»ft) 


■V^r  To;  4*  l  Oo  -4-1  VWU)  ^  v 

i.twl  wn,  if! vii;.  h.|Iu„  4. 


♦YisX 

W't 

.  X. 

(i)-^=^*[cosa  Lw  ^  + 


t-i)] 


(  L 


s,  ’-Jvx  -  %ivX 


A9MOU8  RESEARCH  P  O  U  N  0 A  r  I  O  N  OP  ILLINOIS  INSTITUTE  OP  TECHNOLOGY 

F-14 


W 


umt- 


M 

&, 


'o' 


III, 


'■'J 


+  2UV,<j„+ 

ff  Z^»a)+ 1 1  v;'q„u) 


tf&V  -  #ir-  U-0 

..  l=t  e«  =  <p(iw) 

cfC^  ^  TvlE.  &'.FUkJ$  to.UNT.OM 

‘p'CiKis^V 

<?''(!)*  j|s<PU) 

Taeiu  09=  •.  -■  -• 

<Ml) 


AHMOUH  BESEA9CM  FOUNDATION  Of  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


■  0,  V  04  D-J04.  4  Oc"0^&4, 

^i<fl  =  CQk  ¥&4,Cz) 

An  -  [u',+  Jl-  *[  0A  +•(&,*■  C-s}  ^col  ^  4" 


PROJECTILE  EQUATIONS 


P(.0  2^f  1,  ^  +  *■  ^\i^a  4^1*  ^XQ^uJ'  ^ 

“(A4- 4-2  ^toAta+^^uXt  *-tJL v  ^  L  Li)  * 

(f*i 

4  A  4- ■  A  i*- vT-w  &n  -  04, A(X  4- 

4.  A«o  +^.<e^dfaUa-^»44^u)'] 

Mp^  m  ^p  *  ^CeG^g,*-  ^)  +  (Ag  4- 04, 

-  *(*4*  «\l5 - 0fa(A, - (  U U  *  X -*-£,)  '+ 

4^V  V9A^Ai5)v0^5]ut+  ^  (p"X) 


noun  R  S  S  5  *  «  CD  fOUNOlTIQN  Of  ILLINOIS  INSTITUTE  Of  TECHNOLOGY 


B.  DETERMINATION  OF 

D 

a,.4^&.a,4h.  ■*  t>„) 

,-K  +  «*p)  i&(Alo)  Au.) 

c.  TUBE  EQUATIONS 

^  *lxl  -  *.]  <?« = M) 


Coc^^j-Vyv  )  “C^X-Qc^As^Uj 


-  S»sr 

V,4 

■>  f  ^  4-1 W^4+ fc*  As}  + 

A^.-  ^  (A^Aj,)  +  Be,  A  7^  “^”4."+  (p-fl) 

Y)  Pv  f  ^  ^  ' 

“  %“P5  ^4-9toUo-A5)+6^x](u+l7)  + 

+”(^3  +  OfeAx^J-Afe'C^fc  4  <j)  4  &.J  | 

^  +[C  (a6  -  2  ofa  Atfa + 4a;)}  Vi  n» 

(’i)-Cy'\r,'^)]W(|)\4  (L-^  -v- 
4  Uv  &<•  &s)[  ^  ~  ^  Cos.^3  ^ 


(p-9) 


a[  A,^- i\'Q 


k+.~.+g:  U  + 


'!  K 


>  t-v>r 

b4  tAvv  ^  An  4  (&cAb"9 


-  U«,  ^  19*.)  ill 


4  A  Is,,  4 


*»MOUR  »£S»»C«  P  O  L)  N  0  A  T ,  O  N  OP  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


D.  END  SHEARS  AND  MOMENTS 


1 

MitOjt')  »-£M.cos.*v  Ell 
1>*M  =  Eil5*k, 


My(ci,0=  -til  A„C^  J 

E.  BREECH  EQUATIONS 

C-fL  V  (A4  *"/) 

v  4  An]  = 

=  fLCti 1  Ll)^4.  -p**  -fu) 

1  4^4.4^)  ^  + 

A  (Ac*-^)  101  a 
-  “?Cy-S^(0^)-[  l-»Ul)3KJw 

W'b[-9a  (vjft  ^  -  A,-,  a  -  Ag  jtf,  4  Au.  V  (  Atto-  ^0) 

v  u^  aOI^c 

“  ^c.?r 

ARMOUR  9ESEARCH  fOUNO*flON  OF  ILLINOIS  INSTITUTE 


(F- 10^ 


(p-u) 


lP-a) 


(F-13) 


OF  TECHNOLOGY 


G.  TRAILS  AND  GROUND  REACTIONS 


t-IM  IT  *  O 


e  y-°  u-16) 

Ton-  a-  >o 


*■  ♦‘fc')  Foe.  ^70 

Ift  *  lt»0u 


(F-19) 

Cf-Xo) 

(F-a) 


H.  VERTICAL  EQUILIBRIUM  AT  BALL  JOINT 


£.*V  Aij  2<**  4-  Ac,2wy  ^ 

+•  (&,  +  (k)k<.*  4  t-Gy  *  O 


Cf-ii) 


ANGULAR  MOTION  OF  THE  BALL  JOINT 

~^5  ©4-  *■  jit.  ->‘w  4-  Tu*(25  -Tto^Ss,  V: 

^MOe^o 

"1(^1  ^4  "  ^4-  63-  t^r.5.  6,  4-Ti,*  o  s  Ls  ^ 

-(i>c a^)  4*  -  o 


(y-z^ 


V-a0 


unoue  BSSTAecH  touhoition  of  iuinois  instituts  of  hchnoloot 


F -22 


VII.  COMPILATION  OF  UNKNOWNS 


Unknown 

Equation 

Unknown 

Equation 

Unknown 

Equation 

E6x 

631 

yT 

619 

yB 

621 

3 

623 

5 

627 

t 

T6y 

626 

R6x 

628 

R6* 

630 

E6y 

632 

Pcy 

611 

N 

cx 

613 

N 

cz 

615 

607 

6 

606 

C 

y 

604 

600 

N 

P* 

602 

2 

620 

1 

624 

4 

622 

T6x 

625 

P 

cx 

616 

R, 

6y 

629 

u 

610 

P 

cz 

612 

N 

cy 

614 

S  ,  S 
y  * 

609 

M  ,  M 
z  y 

h 

n 

608 

CP« 

603 

C 

z 

605 

N 

py 

601 

A8MOU8  SE5EKCH  FOUNDATION  OF  ILLINOIS  INSTITUTE  OF  TECHNOLOGY 


DIAGRAM 


APPENDIX  G 
FIXTURE  DESIGNS 

I.  INTRODUCTION 

The  purpose  and  some  of  the  operations  of  the  experimental  structure 
modifications  are  described  in  the  body  of  the  report.  This  appendix  shows 
the  details  of  the  construction  of  the  more  important  fixtures,  the  dampers, 
and  the  elevating  system  flexibilities.  Not  only  the  springs  and  dampers 
used  on  the  program  are  shown,  but  also  a  design  showing  a  combination 
spring  and  damper  which  could  be  used  in  the  actual  launcher.  (The  fixture 
used  in  the  experimental  program  prohibited  changes  in  elevation.) 

II.  FIXTURE  DESCRIPTION 

Figure  G>1  shows  the  assembly  of  the  auxiliary  dampers  used  for  the 
firing  program.  It  is  simply  a  piston  submerged  in  a  fluid  inside  a  cylinder; 
the  piston  has  orifices  so  th»t  its  motion  gives  a  force  approximately 
proportional  to  the  square  of  the  velocity.  A  detail  drawing  of  the  piston  is 
shown  in  Figure  G-2. 

Figure  G-3  shows  the  assembly  of  the  flexibility  added  to  the  elevating 
rod;  Figure  G-4  gives  the  spring  characteristics.  This  assembly  was  inserted 
into  the  present  elevating  system  simply  by  moving  the  existing  elevating 
screw  up,  and  placing  the  spring  assembly  between  the  support  in  the  upper 
carriage  aide  and  the  screw.  The  flexibility  could  be  removed,  if  desired, 
by  replacing  the  spring  by  a  rigid  spacer.  An  alternate  method  was  used, 
however,  that  of  replacing  the  entire  elevating  screw  by  more  rigid  rods. 

Figure  G-5  shows  a  possible  configuration  of  a  spring-damper 
combination.  The  entire  spring  is  submerged  in  hydraulic  fluid;  in  this  case 
the  spring  end  retainers  have  orifice  holes  and  act  as  the  damper  pistons. 

The  lower,  smaller  spring  and  piston  assembly  is  a  device  which  acts  as  a 
fluid  reserve  which  compensates  for  changes  in  ambient  temperature.  The 
entire  assembly  would  be  keyed  so  that  it  could  transmit  the  torques  for 
elevating  the  launcher. 
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MATL.  •  •  •  . . AISI  1095  SPRING  STL. 

WIRE  DIA.  . . . 1/2  " 

MEAN  DIA. . .  1-7/16" 

ACTIVE  COILS . 4-1/3 

TOTAL  COILS  •••’•••••  6-1/3 

RATE  .  7000  lb/ In. 

FREE  HEIGHT  •  •  •  • . 3-13/16 

SOLID  HEIGHT .  3-3/16 

INSTALLED  HEIGHT . 3-5/8 

PRE-LOAD . 1  ZOO  lb 

DEFLECTION . .572"  '  .j  . 

MAX.  LOAD  •  ■  . . .  •  4000  lb 

STRESS . .  115,  000  p.  s.  i. 

WOUND  •  •  . . R.H.  -  L,  H. 

ENDS  •  •  . . .  SQUARED  GROUND 


Fig.  G-4  SPRING  DATA 
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APPENDIX  H 

STRUCTURAL  STIFFNESS  MEASUREMENTS 
R.  M.  Drach  and  R.  H,  Van  Beek 

I.  INTRODUCTION 

In  order  to  improve  launcher  accuracy  by  choosing  optimum 
launcher  parameters,  the  individual  effect  of  each  significant  parameter 
must  be  determined.  Studies  with  mathematical  models  are  being  used  to 
determine  these  effects,  followed  by  firing  programs  to  substantiate  the 
theoretical  conclusions.  The  success  of  using  mathematical  models  for 
this  purpose  depends  upon  at  least  two  major  points:  (1)  the  degree  to 
which  the  models  qualitatively  represent  the  dynamics  of  the  launcher  and 
(2)  the  quantitative  agreement  between  the  mathematical  model  output  values 
and  actual  measurements  of  launcher  motion.  In  order  to  optimize  the 
quantitative  agreement,  the  numerical  input  values  to  the  mathematical 
models  should  correspond  closely  to  measured  values  for  the  actual  launcher. 
For  this  reason,  static  stiffness  measurements  were  made  of  a  number  of 
launcher  structural  components  using  Prototype  No.  3  as  well  as  the  trails 
of  Prototype  No.  1.  Another  equally  important  reason  for  making  these 
measurements  is  to  enable  comparison  of  the  actual  stiffness  values  of  the 
manufactured  parts  with  the  intended  design  stiffness. 

The  overall,  or  combined,  stiffness  of  the  launcher  structure 
below  the  trunnions  is  desired.  It  would  also  be  very  useful  to  determine 
which  components,  such  as  the  carriage  box-sections,  trail  hinge  -pins.  etc.  , 
contribute  the  most  to  the  total  flexibility.  As  it  would  have  been  time  con¬ 
suming  and  costly  to  disassemble  the  entire  structure  and  measure  the  stiff¬ 
ness  of  each  part,  the  stiffness  of  the  entire  structure  was  first  measured, 
and  then,  wherever  disassembly  was  simple,  as  for  example,  the  removal 
of  the  trails,  separate  measurements  were  taken.  An  attempt  was  made  to 
determine  the  stiffness  of  certain  components  by  eliminating  their  contri¬ 
butions  to  flexibility  by  means  of  auxiliary  stiffeners.  In  general,  no 
special  consideration  was  given  to  the  measurement  procedures  to  develop 
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extreme  accuracy  since  stiffness  values  accurate  within  5%  to  10%  were 
regarded  as  sufficient  with  respect  to  the  mathematical  models.  Standard 
laboratory  test  machines  were  used  wherever  possible. 

IX.  MEASUREMENTS 

Static  stiffness  measurements  were  made  for  the  following 
components  and  assemblies: 

1.  The  entire  structure  below  the  t'runnions( supporting  structure). 

2.  The  entire  structure  below  the  trunnions  with  the  upper  and 
lower  carriage  box -flections  clamped  together. 

3.  The  entire  structure  below  the  trunnions  with  the  carriuge 
box-sections  clamped  together  and  struts  inserted  between 
the  trunnions  and  trails.  (Trails  were  parallel  to  the  firing 
tube.) 

4.  The  entire  structure  without  stiffeners>but  traversed: 

a.  20°  right 

b.  20°  left 

5.  Elevating  screw  support 

6.  Trails: 

.  a.  Prototype  No.  3 

b.  Prototype  No.  1 

The  stiffness  of  the  entire  supporting  structure  was  measured 
by  applying  a  known  load  perpendicular  to  a  bar  supported  through  the 
trunnions,  and  measuring  the  deflection  of  the  trunnions  with  dial  gages. 

A  side  elevation  of  this  system  is  shown  in  Fig.  H-l.  The  load  was  meas¬ 
ured  with  calibrated  strain  gages  on  the  loading  bar. 

The  trail  stiffness  was  measured  with  the  two  trails  clamped 
together  by  a  relatively  stiff  fixture  at  their  hinge-pin  ends.  This  enabled 
the  measurement  and  loading  of  the  system  to  be  performed  as  a  simply 
supported  beam  (see  Fig.  H-2b).  Again,  the  deflections  were  measured  with 
dial  gages. 

The  load-deflection  characteristics  of  the  jack  screw  supports 
in  the  upper  carriage  sides  were  measured  by  the  procedure  shown  in 
Fig.  H-2a.  The  compressive  load  applied  to  the  screw  was  determined,  as 
was  the  displacement  of  the  screw  relative  to  the  upper  carriage. 


One  parameter  required  by  the  mathematical  models  is  a 
dynamic  equivalent  for  the  mass  moment  of  inertia  of  the  supporting  structure, 
about  a  line  through  the  ball  joint  perpendicular  to  the  plane  of  elevation.  This 
equivalent  inertia,  which  permits  the  elastic  carriage  to  be  represented  by  a 
single  flexibility  and  single  inertia,  can  be  measured  indirectly  by  measuring 
the  first  mode  natural  frequency  of  the  elastic  carriage.  In  order  to  do  this, 
an  electro -dynamic  shaker  was  attached  to  the  bar  through  the  trunnions  for 
the  purpose  of  applying  to  the  launcher  a  harmonic  force  of  variable  frequency. 
The  dynamic  stiffness  was  measured  as 


where  F  is  the  applied  force,  w  is  the  frequency,  and  a  is  the  acceler¬ 
ation.  Zero  or  very  small  values  of  the  dynamic  stiffness  indicate  reso¬ 
nance,  and  the  smallest  of  the  resonant  frequencies  is  the  first  normal  mode 
frequency. 

HI.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

Figure  H-3  shows  the  load -deflection  characteristics  for  the 
carriage.  When  converted  to  angular  rotation  about  the  base  pivot,  the 
average  of  the  two  aides  yields  a  stiffness  for  the  carriage  of  32  x  10^  lb-in.  / 
radian.  Figure  H-4  shows  the  stiffness  of  the  same  assembly,  but  with 
wedges  between  the  upper  and  lower  carriages  and  links  of  ties  bolting  the 
two  together.  The  stiffness  is  38.  5  x  10^  lb-in.  /radian.  Figure  H-5  shows 
the  stiffness  of  the  same  assembly  as  Figure  H-4,  but  with  additional  links 
connecting  the  trunnions  to  the  trails.  The  stiffness  is  55.5  x  10^  lb-in. /radian. 

Figures  H-6  and  H-7  show  the  load-deflection  characteristics  of 
the  carriage  in  full  traverse,  right  and  left,  respectively.  Traverse  to  the 
right  stiffens  the  right  trunnion  and  adds  flexibility  to  the  left  trunnion,  and 
vice  versa;  but  there  is  a  net  decrease  in  the  average  stiffness.  For  both 
right  and  left  traverse,  the  average  stiffness  is  25  x  10  lb-in.  /radian. 

Figure  H-8  shows  the  load-deflection  characteristic!}  of  the 
lower  elevating  screw  supports  loaded  in  compression.  Thi3  deflection 
occurs  in  the  complicated  structure  consisting  of  the  elevating  screw  gear 
box  and  its  mounting. 
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Figure  H-9  shows  the  load-deflection  characteristics  of  the 
Prototype  No.  3  trails  joined  at  the  hinge-pins  and  loaded  as  a  single,  simply 
supported  beam.  Converted  to  angular  rotation,  which  would  be  produced 
about  the  base  n.'vjt,  the  combined  stiffness  of  the  two  trails  is  145  x  10^ 
lb-in, /radian.  Figure  H-10  gives  the  ;;ame  data  for  the  Prototype  No.  1 
trails,  for  which  the  combined  stiffness  is  109  x  10^.  The  trails  on  the  two 
prototypes  are  of  different  lengths,  and  the  stiffness  varies  as  the  length 
squared.  For  purposes  of  comparison,  the  stiffness  of  the  Prototype  No.  3 
trails  was  extrapolated  to  the  value  they  would  have  if  their  length  were  the 
same  as  that  of  the  Prototype  No,  1  trails.  This  value  is  223  x  10^  lb-in.  /radian 
The  weight  of  the  Prototype  No.  3  trails  is  385  lb;  of  the  Prototype  No.  1  trails, 
351  lb.  On  the  basis  of  these  stiffnesses  and  weights,  the  stiffness -to -weight 
ratio  of  the  Prototype  No.  3  trails  is  computed  to  be  87%  greater  than  that  of 
the  Prototype  No.  1  trails. 

Figure  H-ll  shows  the  results  of  the  dynamic  stiffness  measure¬ 
ments  on  the  carriage.  The  data  are  not  conclusive,  but  it  is  probable  that 
resonance  in  the  first  normal  mode  is  indicated  at  70  cps.  The  inertia  of  a 
torsional  system  having  a  stiffness  of  32  x  10^  lb-in,  /radian  and  a  natural 
frequency  of  70  cps  agrees  reasonably  well  with  the  calculated  equivalent 
inertia  of  the  carriage. 

Each  el^vaHug  *rr»w  and  upper  bracket  assembly  was  loaded 
in  compression  in  a  jniversal  testing  machine.  Fig.  H-12  shows  the 
resulting  load-deflection  curve  for  each  screw-bracket  assembly.  The 
nonlinearities  are  due  primarily  to  nonlinear  local  deflections  under  the  balls 
in  the  ball  nuts  and  also  to  clearances  between  the  balls  and  the  nut. 
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ig.  H -5  STIFFNESS  OF  CARRIAGE  AT  TRUNNIONS  WITH  WEDGES.  TIES, 
AND  LINKS  BETWEEN  TRUNNIONS  AND  TRAILS 


H-6  STIFFNESS  OF  CARRIAGE  AT  FULL  TRAVERSE.  20  RIGHT 


*u;  M4ul}30{J9(l  Iio^uunjj, 
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Fig.  H-8  STIFFNESS  OF  ELEVATING  SCREW  SUPPORTS 


H-15 


NATURAL  FREQUENCY  OF  LAUNCHER  CARRIAGE 


APPENDIX  I 
FIRING  PROGRAM 


This  Appendix  present*  a  list  of  the  firings  conducted  in  order 
to  obtain  the  launcher  response  in  the  vertical  plane  (see  Appendix  J).  The 
instrumentation  for  each  firing  is  listed  as  well  as  the  structural  configura¬ 
tion.  The  bursts  are  coded  with  B  followed  with  a  number;  each  round 
is  individually  marked  by  R  with  a  corresponding  number. 

Because  of  instrumentation  difficulties,  some  bursts  were 
repeated.  These  are  not  listed  because  they  were  under  identical  conditions 
Also  some  additional  firings  were  made  such  as  those  to  study  the  horizon¬ 
tal  dispersion  problem;  these  are  doscribed  in  the  report  body  and  not 
listed  here. 
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Burst  Round  Elevation  [Target  Instrumentation  Fixtures  Movies  Comments 


The  following  fixtures  were  required  for  the  experimental 
firings  of  Pilot  No.  3  at  Ft.  Sheridan.  The  corresponding  numbers  in  the 
following  list  will  be  used  in  the  firing  schedule. 

1.  Variable  flexibility  elevating  rods  with:  a.)  coil  springs  and  b. )  friction 
damped  ring  springs. 

2.  Rigid  elevating  rods, 

3.  Hydraulic  dampers  between  carriage  sides  and  cradle  with  interchange-* 
able  pistons, 

4.  Shorter  pad  placement  on  the  trails, 

5.  Trunnion- to -trail  stiffening  struts, 

6.  Brackets  and  wedges  for  stiffening  the  carriage  box-sections. 

The  following  instrumentation  was  required  for  the  experimental 
firings  of  Pilot  No.  3  at  Ft.  Sheridan.  The  combination  of  instruments  to 
use  for  each  particular  firing  will  be  found  on  the  firing  schedule.  Each  of 
the  following  list  will  be  referred  to  by  number. 

1.  One  strain  gage  on  the  rear  plate  of  the  recoiling  parts  for  a  firing 
signal. 

2.  One  strain  gage  on  each  trail. 

3.  Two  linear  pots  between  cradle  and  ground,  one  on  each  side  of  the  *• 

civkll *i,  perpendicular  to  it  at  45  elevation,  «.  few  inches  behind  the 
cradle  front.  Maximum  expected  motion,  +  3  in. 

4.  One  strain  gage  on  each  elevating  rod  including  the  fixtures.  Maximum 

r  300  x  10  in. /in. 

5.  Two  linear  pots,  one  on  each  side  of  the  launcher,  trunnion  to  ground, 
parallel  to  the  cradle  at  45°  elevation.  Expected  maximum  motion, 

+  3/4  in. 

6.  Two  linear  pots,  one  on  each  side,  trunnion  to  ground,  perpendicular  to 
the  cradle  at  45°  elevation.  Maximum  expected  motion,  +  i/2  in. 

7.  Two  strain  gages,  one  on  each  damper  rod,  to  measure  tension  and 

compression.  Maximum  expected  strain,  -,nn  ,  _-6  -  , ■ 

^  v  200  x  10  in.  /in. 

8.  One  linear  pot  on  each  trail  end.  perpendicular  to  the  ground.  Maximum 
expected  vertical  displacement  (  +  6,  -1/4)  in. 

(The  firing  base  will  be  staked  to  prevent  rearward  displacements.) 

9.  Strain  gages  on  the  indexing  system  shafts. 

10.  Recoil  system  pressure  (one  gage,  right  or  left  side) 

INSTRUMENTATION  NOTE:  Gages  1,  2  and  II  (rear  plate,  trail,  &  recoil 

pressure)  will  be  used  for  all  firings  unless 
otherwise  noted. 
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APPENDIX  J 

EXPERIMENTAL  LAUNCHER  RESPONSE  CURVES 


I.  INTRODUCTION 

The  list  of  the  firings  and  the  instrumentation  for  each  burst 
in  Appendix  I  shows  that  a  very  large  amount  of  data  was  recorded.  The 
most  significant  portion  of  this  data  is  that  which  actually  illustrates  the 
response  of  the  launcher.  In  order  to  limit  the  amount  of  data  to  a  reason¬ 
able  level,  only  the  cradle  response  and  trail  strain  are  shown.  (Exceptions 
to  this  are  Fig.  22  and  23  in  the  body  of  the  report. )  These  curves  are 
presented  for  each  of  the  structural  configurations. 

H.  DISCUSSION 

The  cradle  response  is  shown  in  inches  of  displacement.  This 
is  the  motion  of  the  cradle  front  perpendicular  to  the  plane  containing  the 
trunnions  and  firing  tube.  If  the  trunnions  remain  motionless,  this  cradle 
motion,  by  dividing  by  the  proper  distance,  gives  cradle  rotation.  For  this 
purpose,  the  trunnion  motion  in  the  vertical  plane  was  found  experimentally 
to  be  negligible;  dividing  the  cradle  displacement  by  60  in.  is  a  good 
approximation  of  cradle  rotation. 

If  the  trails  respond  primarily  in  their  first  mode  shape,  the 
trail  strain  is  approximately  proportional  to  the  angular  deflection  of  the 
trail  at  its  connection  with  the  lower  box  member.  Again,  experiments  show 
that  the  use  of  the  trail  strain  to  indicate  rotation  is  a  good  assumption.  This 
strain  and  the  cradle  displacement  serve  to  show  the  launcher  response. 

V  The  following  figures  show  the  cradle  response  and  trail  strains 

: or  the  birat  conditions  described  in  Appendix  I. 
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Fig.  3-2  CRADLE  MOTION.  BURST  B -2 


'utvaic  ti*ji  f 


J  ■*'  .  N  O  —*  N  « 

•ui  ';ujui43«{d«i{j  |u94j  »tp»J3 


SESCADCH  FOUNDATION  Of  ILLINOIS  INSTITUTE 


»•#  »•« 


TRAIL  STRAIN.  BURST  B-7 
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